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ABSTRACT 

We present new observational determinations of the evolution of the 2-10 keV X-ray 
luminosity function (XLF) of AGN. We utilise data from a number of surveys including 
both the 2 Ms Chandra Deep Fields and the AEGIS-X 200 ks survey, enabling accurate 
measurements of the evolution of the faint end of the XLF. We combine direct, hard 
X-ray selection and spectroscopic follow-up or photometric redshift estimates at z < 
1.2 with a rest-frame UV colour pre-selection approach at higher redshifts to avoid 
biases associated with catastrophic failure of the photometric redshifts. Only robust 
optical counterparts to X-ray sources are considered using a likelihood ratio matching 
technique. A Bayesian methodology is developed that considers redshift probability 
distributions, incorporates selection functions for our high redshift samples, and allows 
robust comparison of different evolutionary models. We statistically account for X-ray 
sources without optical counterparts to correct for incompleteness in our samples. 
We also account for Poissonian effects on the X-ray flux estimates and sensitivities 
and thus correct for Eddingtion bias. We find that the XLF retains the same shape 
at all redshifts, but undergoes strong luminosity evolution out to z ~ 1, and an 
overall negative density evolution with increasing redshift, which thus dominates the 
evolution at earlier times. We do not find evidence that a Luminosity-Dependent 
Density Evolution, and the associated flattening of the faint-end slope, is required 
to describe the evolution of the XLF. We find significantly higher space densities of 
low-luminosity, high-redshift AGN than in prior studies, and a smaller shift in the 
peak of the number density to lower redshifts with decreasing luminosity. The total 
luminosity density of AGN peaks at z = 1.2 ±0.1, but there is a mild decline to higher 
redshifts. We find > 50% of black hole growth takes place at z > 1, with around half 
in L x < 10 44 erg s" 1 AGN. 

Key words: galaxies: active - X-rays: galaxies - galaxies: evolution - galaxies: lu- 
minosity function, mass function. 
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Active Galactic Nuclei (AGN) are an important constituent 
of the Universe, playing a vital role in the formation and evo- 
lution of galaxies and having a wider influence on the struc- 
ture of the Universe. Determining the distribution and evo- 
lution of AGN accretion activity throughout the history of 
the Universe, traced by the luminosity function, is essential 
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to constrain models of super-massive black-hole formation 
and growth, the triggering and fuelling of AGN and their 
co-evolution with galaxies. This requires large samples of 
objects spanning a wide range of redshifts and luminosities 
to accurately determine the shape of the luminosity function 
and how this evolves over cosmic time. 

X-ray surveys provide a highly efficient method of se- 
lecting AGN, including unobscured and moderately ob- 
scured sources and low-luminosity sources that may not 
be identified as AGN at optical wavelengths. X-ray emis- 
sion is relatively unaffected by absorption, particularly at 
hard X-ray energies (> 2 keV) for moderate column den- 
sities (Nh < 10 23 cm -2 ), and thus provides a direct probe 
of the accretion activity. Large efforts have been made to 
obtain follow-up observations of X-ray sources in various 
X-ray surveys to determine their redshifts, and thus al- 
low the X-ray luminosity function to be determined (e.g . 
IBarger et aDl2003l : ISzokolv et al 1 l2004l : iTrouille et aHl200a i. 
Such studies of the X-ray luminosity function (XLF) have 
revealed AGN are a strongly evolving population, with the 
overall distribution shifting to lower luminosities between 
z ~ 1 and the present day. This evolution can be de- 
scribe d by a Pure Luminosity Evolution pa rameterisation 
(PLE, iBovle et al] Il993l ; iBarger et all 120051 ) in which the 
XLF retains the same shape (a double power-law with a 
break at a characteristic luminosity, L*), but shifts to lower 
luminosities as redshift decreases. However, recent data in- 
dicate that a more complex, Luminosity Dependent Den- 
sity Evolution (LDDE) parameterisation is required to de- 
scribe the evolution of the XLF at b oth low (z < 1) 

and high (z w 1 - 4) redshifts (e.g. lUeda et alj 12003: 

20081 : lEbrero et ail 



La Franca et al.l 120051: I Silverman et'aif |2_ _ 
20091 : lYencho et al.l T2009 ) . In this scheme the shape of the 



XLF changes with redshift, with the faint-end slope flatten- 
ing as redshift increases. This evolution is also characterised 
by a shift in the peak of the space density of AGN towards 
lower redshifts for lower luminosities. 

However, there are remaining uncertainties as to the 
exact form of the evolution, especially at high redshifts. In 
part, this is due to the difficulties of accurately measur- 
ing the faint end of the XLF, where incomplete and un- 
certain redshift informatio n is a serious issue. Th e Chan- 
dra Deep Fields (C DF-N: lAlexander et~ai1 120031 : CDF-S: 
iGiacconi et al.|[2002h have been key in studies of the XLF, 
providing the deepest available X-ray data and substantial 
programmes of spectroscopic follow-up of the X-ray sources 
to determine redshifts. However, even in these fields a high 
fraction (~ 50%) of the X-ray sources remain spectroscop- 
ically unidentified, mainly because the optical counterparts 
of the faintest X-ray sources are also faint (R > 24), and 
thus beyond the limits of current instrumentation. This can 
severely bias determinations of t he XLF, especially at high 
redshifts (e.g. IBarger et al.ll2005h . 

One general approach to address this issue is to deter- 
mine photometric redshifts for the X-ray sources that are 
too faint for optical spectroscopic follow-up. This signifi- 
cantly improves the completeness of the samples, although 
many X-ray sources which lack detectable optical counter- 
parts remain unidentified. However, there are considerable 
uncertainties in such redshift estimates, particularly for the 
faintest optical counterparts, and a risk of catastrophic fail- 
ures in the redshift determinations, which is particularly 



acute for sources at z ~ 1 — 3. It is also necessary to cor- 
rect for remaining incompleteness in such samples, which 
can prove difficult. Additionally, there is a significant risk of 
associating the incorrect optical counterpart with an X-ray 
source, due to the high surface density of optical sources at 
these faint magnitudes, and thus an X-ray source may be 
assigne d completely incorr ect redshift information. For ex- 
ample, IBarger et al.l (|2003l ) were able to assign optical coun- 
terparts to ~ 85% of the X-ray sources in the CDF-N, but 
noted that ~ 25% of the optical identifications at R — 24—26 
are expected to be spurious. Such false identifications of faint 
X-ray sources could significantly bias determina tions of the 
XLF . To reduce these effects many authors (e. g. lUeda et "all 
l2003l : lBarger et al.ll2005l : ISilverman et al.l2008r ) set relatively 
high X-ray flux limits, thus increasing the spectroscopically 
identified fraction and reducing the numbers of optically- 
faint or unidentified sources, but this also reduces the sensi- 
tivity to the key low-luminosity AGN. The small number of 
fields with sufficiently high spectroscopic completeness also 
means most results are susceptible to cosmic variance. 

An alternat ive approach to t he incomp leteness issue 
was adopted by lAird et all l|2008T l (see also iNandra et al.l 
l2005h to determine the XLF at high redshift, where the ef- 
fects can be most severe. The Lyman-break t echnique (e.g. 
ISteidel fc Hamilton! If993l ; ISteidel et ail |2003l ) was used to 
pre-select a sample of objects at z ~ 3 based on their broad- 
band colours in three filters, and deep X-ray observations 
were used to identify AGN. This technique provides an in- 
complete sample (missing heavily reddened sources in par- 
ticular), but the optical selection functions are well-defined 
and were derived using simulations of the optical data. This 
allowed corrections for the incompleteness to be applied, and 
thus the faint end of the XLF at z ~ 3 could be accurately 
measured. The number density of low-luminosity AGN was 
at least as high as previous estimates, and there was no ev- 
idence for the flattening of the faint-end slope at 2 > I that 
is characteristic of LDDE parameterisations. However, this 
work did not account for the uncertainties in the redshift es- 
timates for sources without spectroscopic confirmation, and 
did not present an evolutionary model. 

In addition to uncertainties in the optical follow-up and 
redshift estimates, a number of issues arising from the X-ray 
observations may also affect the determination of the evo- 
lution of the XLF. The faintest sources in Chandra surveys 
are detected with very small numbers of counts, and thus 
Poissonian effects can be significant. Derived X-ray fluxes 
can thus be highly uncertain, and may be underestimate d 
by up to 50% due to the Eddington bias ( Laird et al. 20091) . 
leadin g to uncertainties in the luminosity. Georgakakis et aD 
(2008) presented a new method of determining the sensitiv- 
ity of Chandra observations, which fully accounts for the 
Poissonian nature of the detection limits, and accounted for 
the flux uncertainties and Eddington bias in the determina- 
tion of the X-ray number counts. However, these techniques 
have not been incorporated in previous studies of the XLF. 

In this paper we present a new study of the evolution of 
the XLF, in which we incorporate new deep X-ray observa- 
tions and develop improved methodologies to address some 
of the issues regarding optical counterparts, completeness, 
photometric redshifts and X-ray fluxes described above. Sec- 
tion [2] describes our data, including our samples of hard 
X-ray selected AGN and the data used to perform high- 
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redshift colour pre-selection of AGN. We adopt a likelihood 
ratio method to ensure we assign secure optical counter- 
parts to the faintest X-ray sources. In section [3] we de- 
velop a Bayesian methodology to determine the evolution 
of the XLF, which accounts for uncertainties in redshift es- 
timates and X-ray fl uxes, incorporates t he im proved sensi- 
tivity calculations of iGeorgakakis et al.l (|200ST ) . corrects for 
remaining incompleteness of the optical identifications and 
allows a robust comparison of different evolutionary mod- 
els. Section 13.31 describes how high-redshift rest-frame UV 
colour-selected samples are in corporated into t his method- 
ology (building on the work of lAird et al.ll200sl ). to provide 
a more robust determination of the evolution at z > 1.2. 
Our results are presented in section [4] and discussed further 
in section [5] 

A flat cosmology with S7a = 0.7 and h = 0.7 is adopted 
throughout. 



2 DATA AND SAMPLES 

To determine the evolution of the X-ray luminosity function 
requires large samples of AGN spanning a wide range of lu- 
minosities and redshifts. For our initial investigation of the 
evolution of the XLF our ap proach is similar to most previ- 



ous s t udies of the XLF (e.g . 
120051 ; ISilverman et"ai1 [2008; 



Ueda et al.ll2003l; iBarger et all 
Ebrero et all 120091 ) in that we 



attempt to assign redshifts to all sources in a purely X- 
ray selected catalogue, utilising photometric redshifts to im- 
prove the completeness. Thus AGN are detected over a wide 
redshift range, which is essential for probing the evolution. 
Our methodology (see section [3} requires probability dis- 
tributions that describe the uncertainties in the photomet- 
ric redshift estimates. We therefore adopt the best avail- 
able photometric redshifts with suitable probability distri- 
butions, or generate our own from the available data. To 
obtain the best constraints on the XLF a high fraction of 
spectroscopic identifications is also required. Additionally, 
we require hard X-ray coverage to reduce biases against 
moderately absorbed sources. Three deep Chandra surveys 
have the necessary spectroscopic follow-up, combined with 
the deep multiwavelength imaging required for photometric 
redshifts: the two Chandra Deep Fields and the AEGIS- 
X survey. We supplement this sample with two larger area 
A SCA surveys with near complete spectroscopic follow-up, 
which are needed to constrain the bright end of the XLF. Ta- 
ble [T] summarises the numbers of sources and fraction with 
optical counterparts and spectroscopic redshifts. Further de- 
tails and a description of the photometric redshift estimates 
are given in section 12.11 below. 

For the high-redshift colour preselection approach the 
requirements are different. Deep X-ray data are again re- 
quired, but the soft (0.5-2 keV) band, where the sensitiv- 
ity of Chandra is highest, may be used for selection as this 
probes relatively hard rest-frame energies at z > 2, and thus 
should not be severely biased by absorption effects. Deep op- 
tical imaging data in U n GlZ, or a comparable filter set are 
required for the colour pre-selection, with deep data at U- 
band wavelengths being the major requirement. A descrip- 
tion of the data for the high-redshift work is given in section 



2.1 Hard X-ray selected samples of AGN 

2.1.1 Chandra Deep Field-South 

The original ~ f Ms observation of t he Chandra Deep Field 
South (CDF-S, [G iacconi et al.ll2002T ) has recently been sup- 
plemented with additional Director's Di scretionary Tim e 
to provide a total of ~ 2 Ms exposure l)Luo et alJ 120081 ). 
which we use in this paper. Data reduction and source de- 
tection was carrie d out using the procedures described in 
lLaird et al.l l|2009l ). For our investigation of the XLF we use 
the hard (2-7 keV) selected source catalogue of objects de- 
tected with false Poisson probability < 4 x 10~ 6 , which pro- 
vides a sample of 248 X-ray sources. The X-ray sensitivity 
(area curve) is det e rmine d using the procedure described in 
IGeorgakakis et all l|200Sl ). 

Deep optical data have been obtained in several bands 
in the CDF-S using both ground-based observatories and 
the Hubble Space Telescope. We have sea rched for optica l 
counterparts using the COMBO-17 survey (jWolf et al.ll2004f ) 
source catalogue, in which sources are detected in deep R- 
band images (10cr limiting AB magnitude ~ 25.6) obtained 
with the WFI instrument on the ESO 2.2m telescope at 
La Silla. Counterparts were foun d using the implementa- 
tion of the likelihood rati o method ((Sutherland fc Saunders I 
1 19921 ; ICiliegi et all 120031 ) described in lLaird et all l|2009h . 
This method allows us to assign secure optical counterparts 
to the X-ray source, accounting for the X-ray and optical 
positional uncertainties, the a priori probability of a coun- 
terpart with a given magnitude, and the probability of a 
spurious optical counterpart. We restrict our sample to ob- 
jects with likelihood ratios > 0.5, which provides optical 
counterparts for 169 (~ 68%) of the hard-band selected X- 
ray sources. 

A number of spectroscopic pr ogrammes have been 
carried out in the CDF-S. Following iLuo et alJ il200Sl). we 
have co mbined spectros co pic su rv eys from | Szokolv et al.l 



( 20041). iLe Fevre et all d2004l), iMignoli et al 



3koly et al. 
Ll d2QQ5|), 



Ravikumar et all (I2007T ) , iPopesso et al.l (|2008l ) and 



Vanzella et al.l (|2008l ). This provides spectroscopic red- 



shifts for 110 sources (65% of the optical counterparts, 44% 
of the entire hard X-ray selected sample). 

A variety of photometric redshift estimates have been 
published in the CDF- S area, focussed on both the X-ray 
source population (e.g. IZheng et al.l I2004T ) and larger pop- 
ulations of optically s elected sources (e.g. IWolf etalll2004l : 
iMobasher et alj|2004 ). However, we have chosen to use our 
own photometric redshifts es timates using t he Bayesian pho- 
tometric redshift code (BPZ. lBem'teal2000l ). which allows us 
to use the full redshift probability distribution output by the 
BPZ code, and thus account for the significant uncertain- 
ties in photo- z estimates in our determination of the XLF 
(see section [3] below). Such probability distributions are not 
available for the published catalogues. 

Our photo-z estimates were determined using the 
full 17-band photometry for the X-ray sources (with se- 
cure optical coun terparts) from the COMBO-17 catalogues 
jWolf et al.ll2004f ). The BPZ code uses a \ 2 minimisation 
template fitting approach. A set of 6 template spectra 
are used, cons is ting of the four standard templates from 
IColeman et"ail <|l980l ) (E/S0, Sbc, S cd, and Irr) and tw o 
star-bursting galaxy templates from iKinnev et alJ j|l996h . 
Additional points in colour space are calculated by interpo- 
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Table 1. Details of the hard X-ray selected samples/fields. 



Field RA Dec X-ray exposure Survey area" Nx b Fopt c F Bpcc d 

(J2000) (J2000) (ks) (arcmin 2 ) 



CDF-S 


03:32:28 


-27:48:30 


1933.4 


436.2 


248 


68% 


44% 


CDF-N 


12:36:55 


+62:14:18 


1862.9 


436.0 


303 


67% 


55% 


AEGIS-X 


14:17:43 


+52:28:25 


188.7 


1 155.6 


397 


86% 


40% 


ALSS 


13:14:00 


+31:30:00 




20 880.0 


34 


97% 


97% 


AMSS 








294 386.0 


109 


98% 


98% 



"Area covered by both X-ray data and required optical imaging. 
'Total number of hard X-ray selected sources in the area covered. 
c Fraction of hard X-ray sources with secure optical counterparts. 

^Fraction of hard X-ray sources with secure optical counterparts and spectroscopic redshifts. 



lating between the templates. The BPZ code also adopts a 
Bayesian approach, applying priors based on the expected 
fractions and redshift distributions for different templates, 
and marginalising over the templates to produce the final 
redshift probability distribution, p(z). These distributions 
describe the uncertainties in the redshift estimate due to 
both photometric errors and colour, redshift and template 
degeneracies, which often result in multi-modal probabil- 
ity distributions. We note that BPZ does not include AGN 
or QSO templates. Broad-line QSOs are likely to be opti- 
cally bright, and may be more easily identified spectroscopi- 
cally due to their clear spectroscopic features. The remaining 
AGN population is expected to be fainter, with a significant 
fraction of the optical light coming from the host galaxy, 
thus improving the match to the galaxy templates. Uncer- 
tainties due to the mismatch between the observed AGN 
SEDs and the galaxy templates may be partially reflected 
in the redshift probability distributions, although the lack 
of a representative template may result in catastrophic fail- 
ures, and an unrepresentative p(z). Additionally, the priors 
are based on observations of the galaxy population in the 
HDF-N, and thus may not truly represent the AGN popu- 
lation. 

In Figure Q] we compare the available spectroscopic 
redshifts with the best photo-z estimates, given by the 
mode of the p(z) distribution. Error bars on the photo- 
z represent 95% central confidence intervals determined 
from the p(z) distribution, although will not fully repre- 
sent more complicated multi-modal distributions. The es- 
timates generally agree well at z S pec < 1.5, and have a 
Sz = |(z p hot - 2spec)|/(l + Zspcc) ~ 0.18. There are however 
a number of outliers. A high fraction of these are associated 
with broad-line QSOs, in which the optical light is domi- 
nated by the AGN, and thus cannot be accurately described 
with the BPZ galaxy templates. Excluding these sources re- 
duces the error to Sz ps 0.13. The remaining outliers are gen- 
erally at high redshifts (z spC c > 1-2), and are assigned large 
errors, indicating broad p(z) distributions. However, we note 
a systematic bias that assigns z spcc > 1.2 sources to lower 
redshifts. The p(z) distributions are equally affected by such 
bias; accounting for the redshift uncertainty, however broad, 
will not correct for this effect. Indeed, the systematic bias 
indicates a catastrophic failure of the BPZ code to correctly 
fit the observed SEDs with realistic templates at the higher 
redshifts. This failure may severely bias any determination 
of the XLF. We explore this in more detail in section [5~3l 
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Figure 1. Comparison of spectroscopic and template fitting pho- 
tometric redshift estimates in the CDF-S {circles), CDF-N {tri- 
angles) and AEGIS {squares). Error bars are 95% central confi- 
dence intervals derived from the redshift probability distributions. 
At z < 1.0 the photometric and spectroscopic redshifts are gen- 
erally in good agreement. At higher redshifts the photo-z's are 
less reliable, reflected by the larger error bars. A high fraction 
of the outliers are associated with broad-line QSOs {open sym- 
bols), which are poorly fit by the templates. However, there is a 
systematic bias to assign z > 1.2 sources to lower redshifts. 



2.1.2 Chandra Deep Field-North 

The Chandra Deep Field-North (CDF-N) provides a sec- 
ond 2 Ms X-ray field. A point source catalogue based on 
the full 2 Ms e x posur e obtained has been presented by 
lAlexander et al.l (|2003l ). however we use our own reduc- 
tion and source catalogue, again using the procedure de- 
scribed in lLaird et al.l (120091). Our catalo g ue pro bes to lower 
significances than the lAlexander et al.l l|2003h work, and 
allows us to adopt our improved sensitivity calculations 
l|Georgakakis et al1l200cf ) . In this paper we use objects with 
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significant detections (false Poisson probability < 4 x 10 -6 ) 
in the hard band, providing a sample of 303 X-ray sources. 

A vast quantity of multiwavele ngth data are availabl e 
in this thoroughly observed field fe.g. lDickinson et al.ll2003T ). 
To find optical counterparts to our X-ra y sources we use th e 
Hawaii Hubble Deep Field North data ijcapak et al.ll2004h . 
which fully covers the X-ray data with deep observations in 
U, B, V , R, I, z' and HK' bandpasses. Optical counterparts 
were found using the likelihood ratio method, as in the CDF- 
S. This yields secure optical counterparts for 204 (m 67%) 
of our hard-band selected X-ray sources. 

We searched for spectroscopic redshifts for our secure 
optical counterparts in a variety of sp ectroscopic surveys 



that have been carr i ed out in the region llBarger et al.l 



Trouille et ai1l200Sl: IWeiner et aD 120061 : ISteidel et all 



2008 



2003 



Reddv et al.l I2006T ) . This provided spectroscopic redshifts 



166 sources (81% of the X-ray sources with secure optical 
counterparts, 55% of the total hard X-ray selected sample) . 
The spectroscopic completeness is much higher (~ 80%) for 
the brightest X-ray sources (fx > 10~ 14 erg s -1 cm -2 , see 
iTrouille et aill2008l ) . 

As for the CDF-S, we have recal culated photometric 
redshift estimates with BPZ, using the lCapak et alj (|2004l ) 
U, B, V, R, I, z' and HK' photometry. A comparison of 
photometric and spectroscopic redshifts is shown in Figure 
[1] Again, a large number of the outliers are associated with 
broad- line AGN, in addition to high redshift sources which 
are assigned large errors. Excluding QSOs, the accuracy of 
our photo-z's is Sz ps 0.11. A systematic bias at z spcc > 1.2, 
assigning sources to lower redshifts, is present, as found for 
the CDF-S. 



2.1.3 AEGIS 



AEGIS (|Davis et al.l [2007) is an international multi- 
wavelength survey of the Extended Groth Strip, covering 
a ~0.25 by 2 degree strip of the sky. Deep (~ 200 ks) X-ray 
data have been ob tained along the entire length of the strip. 
lLaird et al. (2009) described the reduction of these data, and 
presented a catalogue of 1325 unique point sources, merged 
over all bands. In this paper we restrict our analysis to hard 
band detected sources (false Poisson probability < 4 x 10 -6 ) 
and those which also fall within the the area with very deep 
u*g'r'i'z optical data from the Canada-France Hawaii Tele- 
scope Legacy SurvejQ- This reduces the sample to 397 X- 
ray sources, of which 343 (86%) have secure i! band optical 
counterparts (likelihood ratio > 0.5). 

Significant spectroscopic follow-up is available in the 
AEGIS field, incl uding redshifts from the DEEP2 survey 
l|Davis et al.ll2007h and ad ditional targeted follow-up of X- 
ray sources with MMT (|Coil et all 120091 '). This provides 
spectroscopic redshifts for 157 sources (40%). 

In this field we have chosen to use photometric red- 
shifts from 2 different sources. Huang et al. (in prepara- 
tion) us e an Artificial Neural Net work/training set approach 
fANNz. lCollister fc Lahavll2004r i. This approach effectively 
fits an arbitrary function to the redshift- colour relation, us- 
ing a set of spectroscopic redshifts to train the neural net- 
work. As such the redshifts are only expected to be valid for 



http://www.cfht.hawaii.edu/Science/CFHTLS 



sources with the same range of properties as the training set. 
The H uang et al. work use s an 3.6/xm IRAC selected cata- 
logue (|Barmbv et al.l [20081 ). further subjected to a number 
of colour cuts, and utilises photometry from MMT (Mega- 
cam) as well as IRAC, MIPS, and DEEP2. The neural net- 
work was then trained using the DEEP2 spectroscopic data, 
and photometric redshifts estimated for all sources which 
satisfy the selection criteria. This provides highly accurate 
photo-z's (Sz ~ 0.04), albeit for a restricted sample. ANNz 
redshifts are adopted for 38 X-ray sources in our sample 
(without spectroscopic redshifts). To account for the error 
we must determine a probability distribution for each red- 
shift. For the ANNz estimates we have based this on the 
distribution of spectroscopic redshifts, for a given ANNz 
estimate, using all objects in the Huang et al. work with 
spectroscopic identification, as demonstrated in Figure [2] 

For the remaining sources in our sample with se- 
cure optical counterpar ts we adopt photometric redshifts 
from lllbert et al.l (|2006h . which are based on the CFHTLS 
u*g'r'i'z data. The redshifts are based on template fitting 
and a Bayesian approach, and thus use the same method 
as the BPZ code. The priors are calibrated using spec- 
troscopic identifications of galaxies in the Dl field of the 
CFHTLS. The estimates are less reliable than the ANNz re- 
sults (Sz ~ 0.16), but do provide redshift information for the 
remaining optically identified X-ray sources. However there 
is a systematic bias for sources with z spe c > 1.2 to be placed 
at lower redshifts, as was found for the CDF-S and CDF-N. 
The authors have made available the full p(z) redshift prob- 
ability distributions for every source, which we have adopted 
to account for the uncertainties, though once again will not 
account for the systematic bias. 

2.1.4 ASCA Large Sky Survey 

The ASCA Large Sky Survey (ALSS. IUeda et al.lll999l ) cov- 
ers a contiguous area of ~ 5 deg 2 near the north Galactic 
pole, and contains 34 sources detected in the hard (2-7 keV) 



band with the SIS detector, down to a flux limit 



10" 



erg s _1 cm -2 (2-10 keV). Op t ical i dentifications have been 
presented by lAkivama et alj l|2000t) . 2 of these sources are 
optically identified as galaxy clusters, 1 is a star and 30 are 
associated with AGN, all of which have spectroscopic red- 
shifts. 1 source remains u nidentified. An a rea cu rve for this 
survey was presented by lAkivama et all |2003), which we 
adopt for our work. We note that this area curve does not 
fully account for the Poissonian nature of the X-ray source 
detection, but as a relatively high flux limit was set for the 
catalogue and the X-ray fluxes will be accurately measured 
this is unlikely to have an adverse affect on our bright end 
determination. 

2.1.5 ASCA Medium Sensitivity Survey 

The largest area survey in our study is the A SCA Medium 
Sensitivity Survey (AMSS, lUeda et al.ll200ll ). We include 
sources from the AMSSn sub-sample, selected in the hard 
(2- 10 keV) band, with optical identifications presented 
by lAkivama et all ll2003l L The sample includes 87 X-ray 
sources, 76 of which are AGN, to a flux limit of ~ 3 x 10 -13 
erg s _1 cm -2 , and is 100% optically identified with spec- 
troscopic redshifts available for all the AGN. We include 
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tive colour pre-selection approach f or identifying AG N at 
high-redshift, building on the work of Aird et al. (120081 ). We 
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Figure 2. (Top): Comparison of ANN.Z and spectroscopic red- 
shifts in AEGIS, shown for all sources (red circles) and X-ray 
sources (black squares). The ANN2 redshifts generally provide ex- 
tremely good estimates, although are available for a limited sub- 
set of objects, and are restricted to 2 < 1.5 where there are suffi- 
cient spectroscopic training data. The probability distribution for 
an ANN2 estimate is determined from the distribution of spectro- 
scopic redshifts for a small slice of photometric redshifts (e.g. blue 
dashed lines). (Bottom): The distribution of (z — z p hot)/(l + z ) 
for the ANN2 redshift range indicated in the top figure. The true 
distribution (black histogram) is smoothed (red curve) to provide 
an estimate of probability distribution of true redshifts associated 
with a given ANNz estimate. 



additional sources with redshift identifications from the 
AMSSs sub-sample (Ueda & Akiyama, private communica- 
tion), which includes 20 AGN; 2 sources in this sample re- 
main unidentified. The area curve for t he combined AMS Sn 
and AMSSs samples was presented bv lUeda et all (|2003h . 



2.2 High-redshift sample 

Due to the catastrophic failure of our photometric redshift 
estimates at z > 1.2, we have also investigated an alterna- 



include four of the fields used in I Aird et al.1 (|2008r i: HDF 



N (fully encompassed by the CDF-N X-ray field described 
above), LALA CETUS, Lynx and EGS1. We additionally 
include the large area covered by both the AEGIS-X survey 
and the CFHTLS field D3 optic al imaging. This f ully encom- 
passes the GWS field used in lAird et all (|200St ). and thus 
supersedes those data. Our high-redshift data overlaps with 
our hard X-ray selected samples described in section [2~T1 in 
sect ion [<01 we explain how data from our hard X-ray surveys 
at z < 1.2 and our high-redshift colour pre-selected samples 
are combined to determine the evolution of the XLF. 



2.2.1 Optical data 

The reduction of the optical data for HDF-N, LALA CE- 
TUS, Lynx and EGS1 and the source detectio n and pho- 
tomet ry wi th SExtracto r has b een described in lAird et al.l 
l|200ctl and iReddv et alj ^OOfj ). and we refer the reader to 
those papers for further details. For the CFHTLS we have 
used the publicly released images and catalogues available 
from Canada Astronomical Data Centei0, March 2008 re- 
lease. The catalogues were generated using SEx tractor, al- 
thoug h with a slightly different configuration to lAird et al.l 
(2008). Initial detection was performed in the i'-band (which 
is the deepest band) , an d mag n itude s were determined in the 
other bands using the iKronl l|l98Cf ) style aperture (deter- 
mined from the i' band data). The data were not smoothed 
to match the seeing in different bands. These differences 
will affect the detection of objects and measurements of 
their colours, although given the depth and quality of the 
CFHTLS data is unlikely to adversely affect the colour se- 
lection of candidates. Nevertheless, any errors, biases and 
scatter due to the source detection and photometry proce- 
dure will be accounted for by our simulations to derive the 
optical selection functions (section |3.3.1[) . 

Samples of z ~ 3 Lyman-Break Galaxies (LBGs) were 
extracted from the o ptical catalogues us ing the standard 
selection criteria from lSteidel et all (|2003h : 

G-ft<1.2 
U n -G > G-% + 1.0 

19.0 < TZ < 25.5 (1) 

with the faint limit of TZ^ 25.0 for EGS1. An additional 
sample of objects at sligh tly lower redshifts is sel ected using 
the BX selection criteria IjAdelberger et al.ll2004l )L 



G 


-TZ 




-0.2 


u n 


-G 


> 


G - ft + 0.2 


G 


-TZ 




0.2(U„ — G) + 0.4 


U n 


-G 


< 


G -TZ + 1.0 



(2) 

These objects are "sub-U dropouts" , falling in a thin slice of 
colour space below the Lyman-break criteria, with a selec- 
tion function that peaks around z ~ 2.5 (see section |3. 3. l|l . 
Objects were selected from the CFHTLS optical data 



2 http:/ /www4.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/community/CFHTLS-SG/docs/cfhtls.html 
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using the same criteria. However, there are significant differ- 
ences between the CFHTLS filter set and our U n GTL filters 
used in the other fields. The 7?,-band falls between the r and 
i' filters, thus total 7£-band magnitudes were calculated for 
every object using an approximate colour correction, given 
by 

K = r + 0.05 - 1.08fV - i') (3) 

which was determined from our model colour distributions 
described in section 13.3.11 These estimated ^-magnitudes 
were used to apply the total magnitude limits, consis- 
tent with our other fields. However, for the selection in 
colour-space we used the original u*g'r' magnitudes in place 
of UnGTL, thus avoiding the effects of uncertain colour- 
corrections. The differences in the filters will lead to different 
selection functions, which we determine in section [3.3. II 

2.2.2 X-ray data 

All of our fields have been observed with Chandra, and have 
exposure times > 200 ks. The data have be en reduced us- 
ing our pipeline procedure (|Laird et al.1 2009'). For our high- 
redshift samples we use sources detected in the soft-band 
with significance < 4 x 10~ 6 . Our colour pre-selection limits 
us to sources at z > 2, and at these redshifts the soft-band 
will probe rest-frame energies ~ 2 — 8 keV. Thus our selec- 
tion should be relatively unbiased for the selection of moder- 
ately absorbed AGN (Nh < 10 23 cm -2 ), and is comparable 
with the hard-band selection employ ed at lower redshifts 
(|Cowie et al.ll2003l ; iBarger et al]|2005h . 

2.2.3 Matching the X-ray and optical samples 

We determined which of our colour pre-selected high-redshift 
objects host AGN by matching the optical and X-ray source 
lists. The matching is performed with the entire 1Z- (or i'-) 
selected optical catalogues using the likelihood ratio method. 
The colour selection criteria are then applied to the matches. 
Matching the entire optical catalogue ensures that X-ray 
sources are matched to the most likely optical counterparts, 
and prevents colour pre-selected objects being erroneously 
associated with nearby bright X-ray sources, which have cor- 
respondingly bright, secure opti cal counterparts (cf. the di- 
rect extraction method used in lAird et aUl2008h . 

In table[2]we list the total numbers of X-ray detected ob- 
jects within each of our colour-selected samples. We searched 
for spectroscopic redshifts for our sources as described in 
sections [2T21 and [2T31 for the HDF-N a nd AEGIS-X fields 
respectively, and from lstern et al.l (|2002T ) for the Lynx field. 
Remaining unidentified sources were assigned redshift prob- 
ability distributions, p(z), based on the optical selection 
functions (see section [3.3. 

2.2.4 Low-redshift interlopers 

UnGTL LBG selection at z ~ 3 is known to be ex- 
tremely clean and c ontains very few low redshift interlopers 
jSteidel et al.ll2003t) . The majority of interlopers are Galac- 
tic stars, which are not expected to be bright X-ray sources 
and thus will not contaminate our sample. However, the BX 
selection, which does not rely on very strong features such 
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Figure 3. Distribution of available spectroscopic redshifts for 
the X-ray detected colour-selected objects. Green: BX; red: LBG; 
black: total. The BX sample is contaminated by low-redshift in- 
terlopers, although around half may be identified as low-redshift 
star-forming galaxies and on the basis of their /x//opt ratio (see 
Figure [4] section 



as the Lyman-break to identify candidates, is more suscepti- 
ble to contamination from low-redshift sources, mainly from 
low-redshift (2 < 0.3) star-forming galaxies in which the 
Balmer-break i s confused with the Lyman-q forest absorp - 
tion at z > 2 l|Adelberger et ah I l2004l ; IReddv et al] [2006). 
Indeed, a significant number of low-redshift (z < 1.2) inter- 
lopers are identified within our spectroscopically observed, 
X-ray detected BX sample (see table [2J Figure [3} ■ By exam- 
ining the positions of these sources in the fx — /opt plane 
(see Figure [4| it is clear that 5 of these low-redshift inter- 
lopers occupy a distinct position, associated with optically- 
bright objects (1Z< 23), but with relatively low X-ray fluxes, 
thus log(/x//o P t) < —1.5. These objects are the expected 
bright, low-redshift star-forming galaxies, with X-ray emis- 
sion associated with star-formation processes. We exclude 2 
more X-ray detected BX sources from our analysis which fall 
in the same region and are thus also suspected to be low- 
redshift contaminants. The faint BX source with TZ~ 23.5 
and log (/x //opt) < —1.5, which has a spectroscopic red- 
shift in the correct range, is also excluded as we suspect a 
star-burst origin (see also lAird et al.|[2008l ). 

The 6 remaining interlopers, however, are associated 
with low-redshift AGN, and thus fall within the AGN region 
-1.0 < log(/x// op t) < 4-1.0. While any AGN which have 
been spectroscopically identified as low-redshift interlopers 
are excluded from our analysis, additional contaminants may 
remain within our unidentified BX sample. Based on our 
spectroscopic sample the fraction of low-redshift AGN con- 
taminants may be as high as ~ 25%, although our sample 
is biased in the sense that a higher fraction of bright ob- 
jects have spectroscopic identifications. Indeed, in the large 
sample o f BX-selected galaxi es with spectroscopic redshifts 
used by IReddv et al.l (|2008T l the low redshift contamina- 
tion in the galaxy sample was ~ 59% for bright candidates 
(TZ< 23.5), whereas for fainter candidates was only ~ 9%. 
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Table 2. Summary of the fields used to determine the high-rcdshift XLF using rest-frame UV colour pre-selection. 



Field 
name 


RA 
(J2000) 


Dec 
(J2000) 


Nh" 
(10 20 cm- 2 ) 


X-ray 
expos mc b 
(ks) 


Survey 
area c 
(arcmin 2 ) 


Selection 






NspeJ 


Mow z 9 

z < 1.4 


N h 


HDF-N 


12:36:55 


+62:14:18 


1.5 


1862.9 


149.1 


BX 


1192 


9 


8 


2 


7 














LBG 


292 


3 


2 





3 


AEGIS-X/CFHTLS 


14:17:43 


+52:28:25 


1.2 


188.7 


1155.6 


BX 


17050 


77 


26 


10 


66 














LBG 


6933 


30 


10 





30 


Lynx 


08:48:56 


+44:54:50 


1.0 


186.5 


243.4 


BX 


647 


2 








2 














LBG 


223 


1 


1 





1 


EGS1 


14:22:43 


+53:25:25 


1.2 


177.8 


358.9 


BX 


son 


8 








7 














LBG 


32!) 


7 








7 


LALA CETUS 


02:04:44 


-05:05:34 


2.2 


173.1 


233.3 


BX 


660 


i 








4 














LBG 


114 


3 








3 



"Galactic column density from lDickev &: Lockmanl (l990). 

'X-ray exposure time after good time interval and background are ltering (average of the individual pointings for AEGIS-X). 
c Total survey area covered by X-ray and optical data. 

d Number of candidates from the optical imaging satisfying the colour selection criteria. 
e Number of the colour-selected candidates with secure X-ray counterparts. 
f Number of the X-ray detected candidates with spectroscopic redshifts. 

9 Number of the spectroscopically identified sample which are low redshift interlopers (z < 1.4). 

h Final number of sources in our sample after excluding spectroscopically-identified low-redshift interlopers and additional suspected 
low-redshift interlopers according to their position in the fx — /opt plane (see section \2. 2, 41 Figure [1J. 
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Figure 4. X-ray flux (0.5—2 keV) vs. 7?.-magnitude for the 
BX (green triangles) and LBG (red circles) X-ray detected sam- 
ples. Black circles indicate objects with spectroscopic redshifts, 
2spcc > 1.4; purple diamonds indicate low-redshift interlopers 
with 2 spcc < 1.4. As discussed in section [2.2.41 a number of the 
interlopers are associated with low-redshift star-forming galaxies, 
which are optically bright but very faint in X-rays. We thus ex- 
clude all sources with log(/x//opt ) < —1.5. This will exclude the 
low-redshift interlopers, in addition to bright high-redshift star- 
burst galaxies tha t may also conta minate the AGN XLF at the 
faintest fluxes (see I Aird et"al . 2008). Objects with 2 spec < 1.4 are 
also excluded from our analysis, but interlopers associated with 
low-2 AGN may still contaminate the remaining spectroscopically 
unidentified BX-selected sample. 



Conversely, in the sample of objects spectroscopically iden- 
tified as AGN (including QSOs and narrow-line AGN), the 
contamination fraction was only ~ 12%, and approximately 
consistent over the entire magnitude range. Clearly a much 



larger, near complete sample of spectroscopic redshifts for 
BX/X-ray selected AGN is required to truly determine the 
interloper fraction within our sample, and could also allow 
improvements to the selection criteria (such as introducing 
data from redder wavebands) to be developed to help reduce 
the contamination. However, given the current uncertainties, 
we have chosen not to make further corrections to account 
for the interloper fraction within our sample. We thus note 
that the faint end of our XLF based on the BX sample may 
be over estimated due to low-redshift contamination. 



2.2.5 High luminosity sample 

The total area coverage for our high-redshift sample remains 
relatively small, and thus we require an additional sample 
of high-luminosity AGN to constrain the bright end of the 
XLF. We have chosen t o use bright objects [fx > 5 x 10~ 15 
erg s -1 cm -2 ) from the lHasinger et al.l (|2005l ) compilation 
of surveys (excluding the CDF-N to avoid overlap with our 
own analysis). These samples provide a soft-band selected, 
highly complete, spectroscopically selected sample. We have 
included objects with 1.9 < z ^ 3.5, providing 30 high- 
luminosity AGN in our high-redshift sample. These bright, 
complete, spectroscopic surveys may be combined with our 
colour pre-selected samples by assuming a top-hat selection 
function over the chosen redshift range (see section I3.3.2|) . 



3 A BAYESIAN APPROACH TO DETERMINE 
THE XLF 

Our compilation of X-ray surveys provides large samples of 
AGN over a range of redshifts with well-defined complete- 
ness, required to determine the evolution of the XLF. How- 
ever, the faint sources from the Chandra fields are subject 
to uncertainty in both the redshift determination (for those 
sources with photometric redshifts or in the high-z colour 
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preselected samples) and the X-ray flux measurement (due 
to the small numbers of counts from the faintest sources). 
This has led us to develop a Bayesian methodology to make 
accurate inferences about the form and evolution of the 
XLF, accounting for the uncertainties in individual sources. 



3.1 The likelihood function 

The first step is to determine the lik elihood f unctio n, £, 
for which we have followed the work of lLoredol {2004). Our 
individual sources are effectively Poissonian points, drawn 
from the distribution given by the XLF. Thus, for a given 
model XLF, the expected number of detectable sources, A, 
is given by 



A 



J dlog L x j 



j^dzA(L x ,z)<f>(L x ,z\6) (4) 



where the integrals are taken over the entire range of pos- 
sible values. A(Lx,z) is the area of the survey sensitive to 
a source of a particular luminosity, Lx, and redshift, z, and 
thus a particular flux if we assume a single spectral shape 
(power law with photon index F = 1.9). A(Lx, z) effectively 
describes the probability of a source being detected. The lu- 
minosity function, <j)(Lx,z j 8), is the differential number 
density of sources per unit co-moving volume as a function 
of Lx and z, given an assumed model described by the set of 
parameters, 6. is the differential co-moving volume per 

unit area, as a funct ion of z . 

As described in iLoredd l|2004h . the likelihood function 
may be constructed from the product of the probabilities 
of the observed data from the individual sources, given the 
model, and the probability that no other objects were de- 
tected. For a Poisson process the probability of no objects 
being detected is simply e _A . In the presence of individual 
source uncertainties, the probability of the observed data, di, 
from an individual source i, given the model XLF is found 
by marginalising pidi | Lx, z) over the model distribution of 
Lx and z. Thus, 



C({ck} | 0) 



(5) 



e~ x x 



f f d.V'" 

Y\ dlogLx —dz p(di | L x , z) </>(Lx, 2 



where {di} represents the set of observed data from all 
sources, and the product in the second part of the equa- 
tion is over M detected sources. For the idealised situation 
where the redshift and luminosity are perfectly determined, 
p(di | Lx , z) will be a delta function at the true values, and 
equation [6] will generalise to the standard maximum likeli- 
hood estimator. For our sample p(di | L%,z) is the convo- 
lution of the redshift probability distribution p(z) and the 
Poisson likelihood of obtaining the observed X-ray data, N 
counts, 



C{N\s,b) 



(s + b) 1 



-(s + b) 



(6) 



where s is calculated from the hard band (2-10 keV) flux for 
a source of luminosity Lx and redshift z, assuming T = 1.9 
and Galactic Nh- Thus the Poissonian nature of the counts 
from the X-ray sources, the associated uncertainty in the 
flux and the Eddington bias are naturally accounted for 



in our determination of the XLF. Assuming F = 1.9 ap- 
proximately corrects for absorption effects due to intrinsic 
column densities Nh < 10 23 cm~ 2 at z ~ 1, although we 
have not chosen to make a more sophisticated correction 
on a source by source basis (see further discussion in sec- 
tion [S3. Our app r oach b uilds on the log N — log S work of 
iGeorgakakis et all |200S), extended to include redshift in- 
formation, and allows us to use the improved method of 
calculating the X-ray sensitivity. For sources with spectro- 
scopic redshifts we assume delta functions for p(z). We also 
assume delta functions for the luminosities of sources in both 
the hard and soft X-ray selected, large area, bright samples 
(sections 12X11 [2~T31 and [233]) . 



3.2 Correcting for incompleteness 

The likelihood function above does not account for the in- 
completeness of our samples. In our hard X-ray selected 
samples, between ~ 14% (AEGIS) and ~ 32% (CDF-N, 
CDF-S) of the sources lack secure optical counterparts, and 
thus have neither spectroscopic nor photometric redshifts. 
However, we can derive a completeness correction term to 
account for the incompleteness in the likelihood function. 

Figure[SJa) shows the distribution of the hard X-ray flux 
and optical magnitudes of our (optically identified) sources. 
Sources are scattered over a range of optical to X-ray flux 
ratios, approximately over the range — 1 < log fx //opt < 1 
although a small proportion are at lower ratios and a sig- 
nificant number of this hard X-ray selected sample exhibit 
log/x//o P t > 1. At faint X-ray fluxes, a large fraction will 
be unidentified due to their faint optical magnitudes. In Fig- 
ure[5jb) we show the log fx/ j opt distribution of our sources. 
A Gaussian approximation, taking the sample mean, /j, and 
standard deviation, a, is shown in red, and provides a rea- 
sonable description of this distribution. If we apply a strict 
magnitude limit of R = 25.5 (sources fainter than this 
limit are now deemed unidentified), and assume that all the 
unidentified sources would have optical counterparts fainter 
than this limit (reasonable given the very deep optical data 
we have used to search for counterparts), then we can in- 
troduce a correction to the expected number of sources, A, 
given by equation 2] based on the fraction of the log fx/ f opt 
distribution probed at a given flux. Thus, 

A = 

J dlogLx j * %dz J dR A(Lx,z)4>{Lx,z\e)g(Lx,z,R\ti,,a) 
= J dlogLx j %dzA{Lx,z) 4>{Lx,z \ 6) C(L x ,z\ fi,,a) 

where g(Lx, z,R\ /x, a) describes the R- magnitude distribu- 
tion of sources at a given flux, /x(Lx, z), and is derived from 
the Gaussian distribution of log fx/ f out- g(Lx, z, R | fJ,, u) is 
normalised such that 



f 

j — ( 



di? g(Lx, z, R | fi, a) = 1. 



Thus the factor 



C(L X , z\ fj,,a) 



dRg{Lx,z,R\ n,a) (8) 



is a completeness correction, reducing the expected number 
of detectable sources with low X-ray fluxes. 
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improve our completeness correction the log fx //opt distri- 
bution should be determined simultaneously with the XLF, 
accounting for the incompleteness, using the observed data 
(the R- magnitudes of the detected sources). We therefore 
introduce g(Lx,z,R) into the second term in equation [6] 
modifying the likelihood to 

£({d,} | e,n,a) = e- x x 



\[ J dlogix J%dxp(di\Lx,z)4>(Lx,z\0)g(Lx,z,Ri 



li, a) 



where we neglect any uncertainty in the measurement of 
the i?-magnitude of an individual source, Ri, which is a mi- 
nor simplification. The parameters /j, and a are introduced 
as nuisance parameters, and will be determined along with 
the XLF, integrating over their uncertainties when making 
inferences regarding the XLF parameters, 9. In Figure [5p 
the green dashed line indicates the best-fit (mode) values of 
fi = 0.24 ± 0.01 and a = 0.68 ± 0.01 for our LADE evolu- 
tionary model (section |4.2[) . The peak of the distribution is 
clearly shifted to a higher fx/ f opt, indicating the bias in the 
simple estimates. 

3.3 Incorporating the high-z colour pre-selected 
samples 

In section [4] we present results for the evolution of the XLF 
over the full available redshift range based on our hard X-ray 
selected sample only (using spectroscopic and photometric 
redshifts). However, due to the systematic bias in photo- z 
estimates at z > 1.2, we also wish to derive the evolution us- 
ing our high-redshift colour pre-selected samples. Due to the 
a priori colour selection these samples will suffer from large, 
but well-defined incompleteness. Accounting for this incom- 
pleteness, and combining our soft-band selected high-z sam- 
ple with the hard X-ray selected samples at lower redshifts 
requires a number of further modifications to our method- 
ology. 



Figure 5. (a) Hard (2—10 keV) flux against R magnitude (AB) 
for our entire hard X-ray selected sample (black circles). Black 
dashed lines: log/x//opt =-1.0, 0.0 and +1.0, indicating the typ- 
ical range of AGN. Green lines: mean (dashed) and ±lcr (dotted) 
of the best-fit Gaussian log/x//opt distribution (see section fF^J • 
Blue dotted line: R = 25.5 limit for including optical counter- 
parts, (b) Histogram of the distribution of log/x//opt (black), 
and the Gaussian distribution corresponding to the sample mean 
and standard deviation (red curve). This distribution will be bi- 
ased due to incomplete optical identifications, as well as the dis- 
tribution of X-ray fluxes. We account for these effects when fitting 
our evolutionary models (see sections 13.21 and l4t. and derive the 
log fx //opt distribution (green dashed curve), which is shifted to 
higher fx/ f opt ratios. 

However, the observed distribution (the red curve in 
Figure [SJ) will be biased by the very incompleteness we are 
trying to correct for, as sources with R > 25.5 will be missed. 
The extent of the bias will depend on the flux distribution 
of the X-ray sources (and thus the XLF, and the X-ray 
sensitivity), which determines the number of sources with 
faint X-ray fluxes where the problem arises. Therefore, to 



3.3.1 Optical selection functions 

The first step is to determine the optical selection functions 
for the colour selection criteria. Thes e were determined using 
the simulation method described in lAird et al.l 1120081) (se e 
also lSteidel et al.ll 19991 ; iHunt et alj2004lReddv et al.ll2008l ) . 
Three sets of template spectra were determined for differ- 
ent spectroscopic c lassifications (QSO, NLAGN or GAL: see 
ISteidel et al.l 120051 ). and used to generate model [/ n G7?.and 
u*g'r' colour dist ributions. GAL and NLA GN colours were 
generated using a iBruzual fc Charlotj (|2003h model template 
spectrum for a galaxy with continuous star formation and 
age 1 Gyr. Th i s tem plate spectrum was reddened using the 
ICalzetti et al.l (|2000T ) relation for obscuration by dust, with 
various extinction coefficients drawn from the distributions 
based on the observed range of UV sp ectral slopes for GAL 
and NLAGN sources in LBG s amples l|Adelberger fc Steidell 
l200d ; ISteidel et al.ll2002l . l2003l ). Such modelling is not a full 
physical description for the NLAGN, being based only on a 
template spectrum for star-forming processes, but is suf- 
ficient to simulate the observed range of spectral shapes 
and redder distribution of broadband c olours for NLAGN . 
For QSOs we have followed the work of IHunt et ail (|2004) . 



The evolution of the hard X-ray luminosity function of AGN 11 



The template spectra were generated by varying the contin- 
uum slope and Lya equivalent widths of a composite of 59 
QSO spectra ( corre cted for intergalactic absorption) from 
ISargent et ail dl989l ), based on the Gaussian distributions 



given m 



SSS». These model spectra for all three 
types were then redshifted over the range z = 1 — 4 and in- 
tergalactic absorption was applied by simulating a random 
line-of-sight to each object with intervening Lyman-limit ab- 
so rption systems distribu ted according to the MC-NH mode l 
of lBershadv et all l| 19991) . as described bv lHunt et ai1i|2004l ). 
The final spectra are multiplied by the filter transmission 
curves to produce the model colour distributions. Figure [5] 
compares colours for a single model spectrum (for each clas- 
sification) at different redshifts in both U n GlZ and u*g'r' 
colour space. 

The full colour distributions were used to add artificial 
objects into the optical imaging data. Selection functions are 
then determined as the probability of recovering an object 
with a given 1Z magnitude with either the LBG or BX selec- 
tion as a function of redshift. Objects were rec overed using 
our SExtractor routines (see lAird et "all l2008h in the U„G 
1Z fields. The SExtractor configuration file used to generate 
the released catalogues for CFHTLS is available on- lin<E 
thus we used this configuration for those data. Object s were 
selected in the i'-band image and flexible iKronl || 19801 ) style 
apertures were used to measure the magnitudes in other 
bands, and calculate the output colours. Our selection func- 
tions will thus account for the different source detection and 
photometry technique, and the effects of the different filter 
sets. 

Figure [7] compares our final selection functions (derived 
from the simulations) for the different colour criteria, filter 
sets, and optical classifications. There are clear differences 
between the UnGlZ and u* g'r' selections, which we expect 
given the varying wavelength coverage of the filter sets, and 
results in variation in our selection functions. The u* filter 
extends to longer wavelengths than the U n filter, and as such 
the u* — g' colour becomes sensitive to the Lyman-break at 
higher redshifts than the U n — G colour, as can be seen from 
the tracks in Figure [5] Thus the LBG selection functions are 
shifted to slightly higher redshifts when using u* g'r' rather 
than the U n GlZ filter sets. The same effect leads to BX se- 
lection functions which probe to higher redshifts, although 
the low redshift cut-off is relatively unaffected. Thus the BX 
selection functions are generally broader for u* g'r' selection. 

The red continuum of our NLAGN model generally 
leads to a reduced efficiency of the selection functions. This 
effect is particularly seen for the BX selection function for 
UnGlZ colours. The example NLAGN model shown in Figure 
[5] has colours which are very close to the boundary of the BX 
selection box. The intrinsic scatter we apply to our models 
and the photometric scatter introduced via our simulations 
will both lead to incompleteness of the BX selection func- 
tion, which is clearly seen in Figure[7] (bottom centre panel). 
However, the r' filter is significantly bluer than the 1Z band- 
pass, resulting in g' — r' colours which are generally bluer 
than G— 71, and thus the NLAGN BX selection function is 
less adversely affected by the red continuum for u* g'r' selec- 



3 http:/ /wwwl.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/community/CFHTLS-SG/ docs / cfhtlscat.sex 



tion, and has a higher efficiency (cf. top and bottom central 
panels in Figure [7}. The same effect is seen for the GAL BX 
selection functions, although is much milder as the distri- 
bution of GAL colours is not as red as the distribution of 
NLAGN colours. 

The model tracks for the QSO colours in Figure [5] follow 
complicated routes which have very different shapes for the 
different filter sets. This is due the presence of the strong, 
broad emission lines, which will move in and out of the var- 
ious filters at different redshifts. However, the basic charac- 
teristics necessary for our colour selections remain present: 
the Lyman-break leads to a rapid reddening of the U n — G 
(or u* —g') colour at z ~ 3; a characteristically flat spectrum 
is found at lower redshifts. 



3.3.2 Modifying the likelihood function 

The optical selection functions described above can be used 
to modify the likelihood function derived in section 13.21 to 
allow us to combine the hard X-ray selected samples at low 
redshifts (z < 1.2) and our colour pre-selected high-redshift 
samples. Firstly, A(Lx,z), which effectively gives the prob- 
ability of detecting a source with given Lx and z, must be 
modified to account for the 2 different selection techniques. 
For z < 1.2, A(Lx, z) is calculated from the area curves for 
our hard X-ray surveys, converting the hard (observed 2-10 
keV) X-ray flux to an Lx (rest-frame 2-10 keV) and z assum- 
ing r = 1.9 and Galactic Nh, as before. For z > 1.2, we use 
the area as a function of soft X-ray flux, multiplied by the 
average of the combined BX and LBG selection functions, 
Pf(z) for all magnitudes and spectroscopic classifications for 
a field, scaled to a maximum of 1. Thus, 



Pf(z) 



E 



P f (lZ,z)dlZ 



(10) 



A(L x ,z) = Y,MMLx,z))xP f (z) [z>1.2] (11) 

fields 

where Af(fx.(L%, z)) is the area of each field sensitive to a 
soft (0.5-2 keV) flux for a given 2-10 keV luminosity source 
at 2, assuming T = 1.9 and Galactic Nh, for each field. 
This reduces the probability of objects being detected out- 
side of the range covered by the colour pre-selection, and 
thus provides realistic expected numbers at high redshifts 
in the likelihood evaluation (equation 1 1 Op . As we consider 
th e uncertainty in the X-ray flux we can utilise the method 
of iGeorgakakis et all (|2008l ) to determine the X-ray area 
curves. For the high luminosity sample (section I2.2.5[l . the 
equivalent selection function is simply a top-hat function 
over the entire redshift range (1.9 < z < 3.5), corresponding 
to a complete spectroscopic follow-up of the purely X-ray 
selected high-luminosity sample. 

For each of the individual colour-selected sources we 
must assume a redshift probability distribution, Pi(z) for 
which we adopt the selection function for the appropriate 
field and sample (unless a spectroscopic redshift is available 
in which case we adopt a delta function). Thus the differing 
redshift coverage of the selection criteria for different types 
and filter sets will provide some additional information re- 
garding the evolution, and this effect will be explicitly ac- 
counted for by our study. The X-ray flux uncertainty is also 
accounted for in an identical manner as for the low-redshift 
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Figure 6. Example model colour distributions for our GAL (black circles), NLAGN (green squares) and QSO (red triangles) classifica- 
tions. The same model for each type is shown in U n G1Z (left) and u*g'r' (right) colour-space. The model is redshifted over 1.0 < z < 4.0 
Marks show the colours at Az = 0.1 intervals, with larger symbols indicating 2.5 < z < 3.5. The variation between colours through 
different filter sets leads to differences in the optical selection functions (see Figure [7] section f3-3.lt . 



hard-band selected samples, only converting from the soft- 
band flux to 2-10 keV rest-frame luminosities, to determine 
p(di | Lx,z) for each of our colour selected objects. 

Additional weighting factors, Wi, must also be intro- 
duced, which are added as exponents of the individual source 
likelihoods in equation 1101 This factor serves two purposes. 
Firstly, it reduces the weight assigned to sources in our hard 
X-ray samples according to the fraction of their p(z) distri- 
bution which falls below z < 1.2. 



So' 2 Az Pi{z) 



(12) 



f °°dz Pi (z) ' 

Thus, sources will be assigned a weight of Wi = 1 if they 
have a spectroscopic redshift z spcc < 1.2, or if their Pi(z) 
based on the photo-z estimate falls wholly below z — 1.2. 
Sources with z spe c > 1-2 or Pi(z) which is all above z — 1.2 
will be assigned a weight of (effectively removing them). 
Sources with photo-z estimates that span above and below 
z = 1.2 will be given appropriate fractional weights. 

For the colour pre-selected objects, Wi is the ratio of the 
total co-moving volume sampled based on our z > 1.2 esti- 
mate of the area curve (equation [TT} , and the the effective 
co-moving volume sampled, based on the selection function 
at the object's 1Z magnitude. 



1 it 



f-A(L t ,z)^dz 



(13) 



where Pf (IZi , z) is the selection function for both BX and 
LBG selection for a field for the observed 7?.-magnitude. 
Thus for the high-redshift sample, the Wi factors apply ad- 
ditional weight to sources with fainter optical magnitudes, 
which were thus less likely to be selected by the colour pre- 
selection technique. 



Our colour-selected sources are also incorporated into 
the calculation of the fx/ f opt distribution, and the subse- 
quent completeness correction. The derived fx/ } opt consid- 
ers hard-band X-ray fluxes. To achieve this an equivalent 
hard-band flux is calculated for each Lx,z, whereas Lx for 
a source was originally derived from the observed soft-band 
flux. 7i-band magnitudes are derived from the ^-magnitudes 
by applying a small colour correction (based on our model 
templates described above). 

Sources identified as low-redshift interlopers (section 
I2.2.4|l are excluded from the colour-selected samples, but 
may still be included in the main sample to constrain the 
z < 1.2 XLF for the AEGIS and CDF-N fields where our 
surveys overlap. 34 sources which were assigned photometric 
redshifts in our original hard X-ray selected samples are also 
present in our BX or LBG samples. Of these 17 were assigned 
Zphot < 1-2; the remainder were assigned to the higher red- 
shifts expected for objects in our colour selected samples. For 
all of these sources we do not consider the photo-z: distribu- 
tion, and consider them solely with the colour selected sam- 
ples, assigning the associated redshift distribution (based on 
the optical selection functions) and weighting factors. 



3.4 Bayesian model comparison and parameter 
estimation 

Having determined the likelihood function, our knowledge of 
the XLF can be fully described by the posterior probability 
function, 

p(0, n, a | {dj}, H) = £(W 1 g#fj 1 H) (14) 



P({*} I H) 
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Figure 7. Optical selection functions for our BX (green solid lines) and LBG (red dashed lines) colour selection criteria. The top row 
shows the selection functions for the CFHTLS-D3 optical data and u'g'r'i' filter sets. The bottom row shows selection functions in 
the HDF-N (with optical data of comparable depth) using U n G1Z filters. Each column is for a different spectroscopic classification as 
indicated. Significant differences in efficiency and redshift coverage exist between the selection functions for different optical types and 
filter sets (see discussion in section [3.3. It . 



where we have introduced the symbol H to represent a par- 
ticular model hypothesis for the XLF, and ir(8, fi, a \ H) is 
the prior probability we assign to values for the parameters 
9, for model H. The denominator, p({di} | H), is known as 
the Bayesian evidence, Z, and normalises the posterior. 

There are now two tasks at hand. Firstly, we wish to 
determine the best estimates of the parameters, 9, for our 
model XLF, and assign them an error. Standard Monte- 
Carlo Markov Chain approaches can be used to draw sam- 
ples from the posterior, and thus provide estimates of the 
parameters. Such approaches ignore the normalising factor, 
p({di} | H), as it is independent of the parameters. Our 
second task is to compare different evolutionary models. 
Bayesian model comparison is achieved by calculating the 
ratio of the posterior probabilities for different models (e.g. 
iKass fc Raftervll 1995ft . 

pjHx I {di}) = p({di} | H{) Tr(ffi) = Zi Tr(ffi) 

p{H 2 | {di}) P ({d z } | H 2 ) tt(H 2 ) Z 2 tt(H 2 ) 1 ' 

where the ratio of the prior probabilities of the two different 
models, ir(Hi)/ir(H 2 ), can be taken as unity if we have no 
prior knowledge which indicates we should favour one of the 
models. The evidence can be calculated by integrating the 
likelihood over the prior parameter space, thus, 

Z= I d D 9 C({di} | 9) n{9 | H) (16) 



where D represents the number of dimensions for the in- 
tegration (i.e. the number of parameters, now incorporat- 
ing the nuisance parameters /i and a into the parameter 
set 9). To perform this integr ation we have implemented the 
"nested sampling" algorithm (Skilling 2 0041 ). which provides 
an estimate of the evidence, Z, along with a numerical un- 
certainty. 800 "live" objects are used for the nested sampling 
to ensure that the large, multi-dimensional parameter space 
is well sampled. Estimates of the posterior distribution may 
be recovered as a by-product of the nested sampling pro- 
cedure, and thus used to determine the values of the pa- 
rameters describing the XLF and their errors. We choose 
to report the posterior mode (equivalent to maximum like- 
lihood) values of the parameters as our best estimates, and 
la equivalent (68.3%) Highest Posterior Density (HPD) in- 
tervals, which minimise the size of the confidence interval 
for the given confidence level. 

It is worth noting the difference of this Bayesian model 
comparison approach to frequentist methods. The frequen- 
tist approach to compare two models would be to calculate 
the likelihood ratio, the ratio of t he maximum likelihood va l- 
ues for the two different models (|Nevman fc Pearsonlll928h . 
The Bayesian approach, on the other hand, integrates the 
likelihood over the entire prior parameter space. As such the 
Bayesian evidence incorporates knowledge of all the uncer- 
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tainties in the data and the parameter estimation. This ap- 
proach also has the advantage that it naturally implements 
Occam's Razor, favouring a simple model with a compact 
parameter space, thus less free parameters, over a more com- 
plicated one, unless the latter provides a signific antly better 
description of the data l|Berger fc Jefferva 1992h . 



3.5 Binned estimates 

Our analysis of the evolution of the XLF considers each de- 
tected object individually, and does not require any binning. 
Our results consist of our best-fit parameter estimates and 
their uncertainties, and the relative Bayesian evidences for 
the different evolutionary models presented in section [4] and 
tables [3] and [4] However, for display purposes it is useful 
to produce binned estimates which represent the value and 
uncertainties over certain discreet ranges of luminosity and 
redshift. Our binned estima tes are based on the -/V b s /-/V m di 
method l|Mivaii et al.ll200lh . which multiplies the value of 
the best-fit model evaluated at the centre of the bin by the 
ratio of the number of observed sources, iV bs, to the pre- 
dicted number in the bin, iV m di, where 



iVmdl 



AlogL x 



ulm'U / j^A(Lx,z)4>{L x ,z\e) (17) 



for each bin of size A log Lx and Az. With our methodology 
a source is assigned a distribution of luminosities and red- 
shifts, and thus can be included in multiple bins. Therefore 
the observed number of sources, iVobs, is found by summing 
the fractional contribution of all sources to each ALx — Az 
bin, 



N ) 



E 



/AlogL x dl °g L xJ Az ^p(A | Lx,z) 

fZo d loe ix Jo°° ^P(di\L x ,z) 



(18) 



Errors on the binned values are based on approximate Pois- 
son uncertainties on the effective number of objects in a bin. 
Points are only shown for bins with an effective iV b s ^ 1. 



4-1.1 Pure Luminosity Evolution (PLE) 

At first we investigated models in which the XLF retains the 
same shape at all redshifts, but undergoes overall shifts in 
luminosity. We assume the standard double power-law shape 
for the XLF, 



d$(L x ) 
d log L x 



Lx 



+ 



Lx 

L, 



721 -1 



(19) 

where 71 is the faint-end slope, 72 is the bright-end slope, 
L» is the characteristic break luminosity and K noT m is a 
normalisation factor for the overall density. Pure luminosity 
evolution can then be desc ribed as evol u tion of the char- 
acteristic luminosity, Ltfz). lBarger et al.l l|2005l ) found that 
the low redshift (z < 1.2) XLF could be well described by a 
PLE model, where L*(z) is given by 



log L* (z) = log Lo+p log 



1 + z 



(20) 



However, prev i ous s tudies (e .g. Ueda et"all 120031 ; 

ISilverman etall 120081 ; lAird et al.l I2OO8I ) have indicated 
that the strong evolution of L* at low redshifts may not 
continue above z ~ 1. A number of authors have utilised 
models which extend equation [20] with ad ditional quadratic 
or higher order term s in log(l + z) fe.g. iPage et al.lll997l ; 
ISilverman et alJl2008h . Alternatively. lUeda et al] (|2003l ) de- 
scribed a PLE model which has 2 different slopes below and 
above a critical redshift, z c . This parameterisation allows 
the luminosity evolution to flatten off, cease completely, or 
increase or decline at high redshifts. We have adopted this 
form of luminosity evolution, additionally allowing for a 
smooth transition between the 2 slopes, thus the evolution 
of L» is given by 



logL*(z) = log L " log 



1 



+ 



l + z c 



(21) 



where the z c parameter controls the transition from the 
strong low-2 evolution to the high-z form. 



4 EVOLUTION OF THE XLF 

We have used our Bayesian methodology and the nested 
sampling algorithm to compare three different evolutionary 
models of the XLF, using both our hard X-ray selected sam- 
ples over all redshifts (using photometric and spectroscopic 
redshifts), and the combination of the hard X-ray sample 
at z < 1.2 and the high-redshift colour pre-selected sam- 
ples. For each evolutionary model we calculated the best-fit 
(mode) for the parameters, with errors given by the Highest 
Posterior Density interval, and the relative logarithmic evi- 
dence, A In Z. Differences in In Z > 4.6 indicate very strong 
evidence in fa vour of the model with the higher evidence 
|jeffrevslll96l1 ). The results are given in tables [3] and [4] for 
the hard X-ray sample only and the combined samples re- 
spectively. Constant (or log-constant) priors were assumed 
for all parameters over the ranges indicated in table [3] 

4.1 Evolutionary models 

We compared the following evolutionary models: 



4-1.2 Luminosity and Density Evolution (LADE) 

The work of lAird et all i|2008l ) indicated that while the XLF 
retained a similar shape at high redshift, there was an overall 
decline in the number density at the faint end (Lx < L*), yet 
little or positive evolution at the bright end. This behaviour 
cannot be described by PLE. We have modified the PLE 
form described above to additionally allow for some overall 
decreasing density evolution with redshift, by allowing the 
normalisation constant K noTrn to evolve as 

or in {z)=\ogK + d{l + z) (22) 

thus introducing an additional parameter, d. With this pa- 
rameterisation the shape of the XLF stays fixed at all red- 
shifts, thus our evolution consists of an XLF which shifts 
in both luminosity and overall density. We refer to this 
evolutionary model as Luminosity And Density Evolution 
(LADE). This model is similar to the Indepen dent Luminos- 
ity an d Density Evolution (ILDE) described bv lYencho et al.l 
(2009), but with a different form for the density evolution 
and a more complex luminosity evolution that can flatten 
off at higher redshifts. 
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4-1.3 Luminosity-Dependent Density Evolution (LDDE) 

As discussed in section[T] a number of authors have proposed 
that an evolutionary scheme in which the XLF undergoes 
a density evolution which is dependent on luminosity is re- 
quired to accu rately describe the evolution of both soft X -ray 
selected (e.g. iMivaii et a .Il2000l; lHasinger et al.ll2005h and 
hard X-ray selected (e.g. Ueda et al.ll2003l ; Silverman et alj 
2008) AGN. We hav e adopted the parameterisation given 
bv lUeda et all (|2003h . in which a typical double power-law 
form (e.g. equation 1 19 p of the XLF at z — is modified 
by an evolution term which is a function of luminosity and 
redshift, 



i/>(L x , z) = i>(Lx, 0) e(L x , z) 



(23) 



where e(Lx,z) is a power-law function of z, which changes 
between 2 different forms at a cut-off redshift, z c , which 
depends on the luminosity. Thus 



(l + z) ei [z<z c (L x )] 

" Mrms)" 12,1 



where 



Zc* [L X > L a ] 

z c *(%y [Lx<L a] ■ ^ 

We allow all the parameters to vary in our fitting (although 
apply priors over constant ranges), thus the evidence, Z, will 
fully account for the complexity of this parameterisation. 

4.2 Results 

Table [3] gives the best-fit (posterior mode) parameters for 
our three evolutionary models using the hard X-ray selected 
samples only over all redshifts, and values of the Bayesian 
evidence relative to the evidence for the PLE model, A In Z. 
These values indicate the the LADE model, introducing an 
additional parameter, provides a significantly better fit than 
the PLE model. There is also very strong evidence favouring 
the more complicated LDDE parameterisation over either 
the LADE or PLE models. Figure [8] plots the maximum- 
likelihood LDDE model along with binned estimates for a 
range of redshifts. It immediately apparent that our LDDE 
model predicts a much milder evolution of the faint end 
of the XLF than the model presented by ISilverman et al.l 
(2008), with significantly less flattening of the faint-end 
slope. This is most likely due to our inclusion of objects 
with fainter optical magnitudes and correction for remain- 
ing incompleteness in the samples (see further discussion 
in section 15.11 below. However, the majority of the flatten- 
ing of the faint-end slope we do observe is found at z > 1.2, 
precisely where our photometric redshifts become unreliable 
and prone to catastrophic failure, and thus this measure- 
ment of the faint end of the XLF will be biased. Therefore, 
we are unable to confidently rule out the PLE or LADE 
evolutionary models using the results from direct hard X- 
ray selection only. 

In table|4]we present the results for the three evolution- 
ary models using the colour pre-selected samples at high red- 
shifts in addition to the hard X-ray samples at z < 1.2. The 
colour pre-selection approach restricts the sample to sources 
for which we can reliably estimate the redshift. This should 



therefore provide more robust measurements of the XLF at 
high-redshifts that are not subject to systematic biases due 
to catastrophic failure of photometric redshifts. The sample 
size is smaller and incomplete, but due to the well-defined 
selection functions we are able to correct for this incom- 
pleteness. We find strong evidence favouring the LADE or 
LDDE models over the PLE model. We can thus rule out 
this simplest model, consisting of an XLF which only evolves 
in luminosity without changing shape or moving in overall 
density, and is thus unable to reproduce the lack of evo- 
lution of the bright end of the XLF at high redshifts and 
the observed decline at fainter luminosities. However, there 
is only very weak evidence (A In Z — +1.1) favouring the 
most complex LDDE evolutionary scheme over the simpler 
LADE model in which the XLF retains the same shape at all 
redshifts. We therefore conclude that there is no significant 
evidence that an LDDE parameterisation, and the associ- 
ated flattening of the faint-end slope, is required to describe 
the evolution of the XLF. Our LADE model, in which the 
XLF retains the same shape but shifts in both luminosity 
and density, is also able to adequately describe the observed 
evolutionary behaviour. Our best-fit LADE model indicates 
strong luminosity evolution takes place between z — and 
z — 1, but is consistent with a value of p2 = 0, indicat- 
ing that the evolution of L, flattens off at high redshifts, 
and the evolutionary behaviour becomes dominated by the 
overall negative density evolution. 

Our LADE model is plotted in Figure [9] evaluated at 
the centre of a number of redshift bins, along with binned 
estimates (see section using both the hard X-ray sam- 
ples and the high-z colour pre-selected samples. The LADE 
model is in good agreement with the bright end of the XLF 
determined from the hard X-ray selected samples at z > 1.2 
(open circles, not used to constrain this model), but gener- 
ally predicts a higher number density below i*, reflecting 
the systematic failure of the photometric redshifts in this 
regime. At 2.0 < z < 2.5 our binned estimates from the 
colour pre-selection are systematically below the LADE best 
fit (although agree within ~ la). Indeed, the binned data 
appear to indicate that the overall number density drops 
between z ~ 1 and z ~ 2, but then rises again between 
z ~ 2 and z ~ 3 (none of our evolutionary models could de- 
scribe such behaviour). However, it is more likely that this 
behaviour reflects additional incompleteness in the BX se- 
lection, which has not been accounted for. The issues and 
limitations of our methodology are discussed further in sec- 
tion 15.21 below. There is much better agreement with the 
2.5 < z < 3.5 binned estimates, where the sample is larger 
and the XLF is more robustly determined. 



5 DISCUSSION 

Our results for the evolution of the faint end of the 2-10 
keV XLF should be the most accurate and robust to date. 
Our work utilises both of the Chandra deep fields, including 
the full 2 Ms exposure available for each field. The addi- 
tional large area of deep (200 ks) X-ray data provided by 
the AEGIS survey are also essential for an accurate deter- 
mination of the faint-end slope. Our sensitive point source 
detection techniques allow us to probe to the greatest depths 
in these data, yet our Bayesian techniques allow us to fully 
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Table 3. Best-fit parameters and relative evidence or our different evolutionary models using the hard X-ray selected sample only 



Parameter 


Lower limit a 


Upper limit a 


PLE 


LADE 


LDDE 




-7.0 


-3.0 


-4 95+°-°| 


-4 02+°^ 

u — 0.15 


-5 91+°-" 


logL, (or L ) / erg s _1 


43.0 


46.0 


44.82t ;°J 






71 


-0.1 


1.5 


0.70±™ 3 3 


n cro + 0.04 


o.8oi o ;°| 


72 


1.5 


4.0 


9 cro+0.10 


9 rrr+0.12 


2.36l°;ii 


Pi 


3.0 


10.0 


5 lfi+0-35 
°- 1D -0.35 


D ' DU -0.58 




P2 


-4.0 


3.0 


i no+0.22 

l.OO_ 22 


63+ ' 42 
U.DJ_ 42 




ei 


2.0 


6.0 






4 4 o+0.30 
4 - 48 -0.30 


e2 


-5.0 


0.0 






-2-8518:2 


z c (or z c *) 


0,1 


3.0 


n 8Q+ 08 
U.»9_ . 0g 


75+ 14 




d 


-1.5 


0.0 




-0 36+ 03 




log L a / erg s _1 


44.0 


46.0 






45 24+ - 19 


a 


0.0 


1.0 






15+ - 01 
u - io -o.oi 



AlnZ . . 0.0 ±0.3 +20.3 ±0.3 ±32.9 ± 0.3 



"Upper and lower limits on the constant priors for the parameters. 



Table 4. Best-fit parameters and relative evidence for our evolutionary models using the hard-band selected samples at z < 1.2 and 
constrained at high redshifts using the colour pre-selected samples. 



Parameter PLE LADE LDDE 

logtfnorm (or K ) j Mpc" 3 -S.lOtHt -4.53±g-g? -6.08t°-g« 

logL, (orLo) /ergs" 1 44.96±g-gl 44.77t^« 44.42t -° 4 

70+ 02 n 62+ 02 77+ - 01 

71 u -' u -0.02 u - Oi: -0.02 u '"-0.01 

014+O.I3 „ m +0.11 9 o n +0.12 
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account for uncertainty in the flux measurements at these 
faint limits, and account for the Eddington bias. By incorpo- 
rating photometric redshift estimates we are able to achieve 
high redshift completeness (~ 75%), but we have accounted 
for the uncertainty in such redshift determinations. We have 
corrected for remaining incompleteness in our samples using 
the relation between optical and X-ray fluxes, which we have 
constrained based on our observed samples. At high redshifts 
we use colour pre-selection to constrain the XLF and reduce 
biases due to the catastrophic failures of photo- z estimates. 
We have also used a likelihood ratio matching technique to 
ensure we only include robust optical counterparts to X-ray 
sources. 

5.1 Comparison with previous results 

Our work offers a number of improvements over previous 
studies. No previous work has considered the errors in X-ray 
fluxes and detection threshold, and thus will not account for 
Eddington bias. Additionally, no prior studies utilising pho- 



tometric redshifts (e.g. Barger et al.l 120051 ; ISilverman et al.l 
l200Sl ;lEb rero et al. 2009) have accounted for the known un- 
certainties in such estimates, or the possible effects of catas- 
trophic failures which we believe led to a biased result at 
high redshifts in our work (see Figure [5] section T4.2J) . The 
commonly used photometric r e dshift catalogue for the CDF- 
N, presented bv lBarger et al.l l|200of) . used the same optical 
data as our work, and the BPZ code, and thus these esti- 
mates will be subject to the same uncertaint ies and potential 
for catastrophic failures. IZheng et al] l|2004h used a range of 
photo- z codes, incorporated a wider range of optical filters, 
deep HST imaging data, and used the X-ray data as an in- 
dication of absorption properties, to provide improved pho- 
tometric redshift estimates of X-ray sources in the CDF-S 
(1 Ms data). Although these estimates were found to be ac- 
curate to within ~ 8%, the authors discussed the remaining 
uncertainties, particularly for the fainter sources, and the 
possibility of catastrophic failure. The lack of spectroscopic 
follow-up makes it difficult to assess the reliability of all pho- 




Figure 8. Our best fit Luminosity-Dependent Density Evolution model of the XLF using the hard X-ray selected samples only, evaluated 
at the centres of the displayed redshift bins (solid lines), along with binned estimates (circles, see section fe-St. The LDDE model ev aluated 
at z = is shown in every panel (dashed lines). We compare our results to the LDDE model presented bv lSilverman et al.l l|2008l . dotted 
lines), which exhibits significantly greater flattening of the faint-end slope at high redshifts. However, we also note that our model fit 
will be biased at z ~ 1.2 — 3 where our photometric redshift estimates are prone to catastrophic failure. 



tometric redshifts at high-z, and the uncertainties should be 
accounted for in estimates of the XLF. 

Our completeness corrections allow us to account for 
the lack of redshift information for a fraction of our sources, 
which is particularly important for determining the faint end 
of the XLF, where the optical counterparts are often too 
faint to be detected in even the deepest available imaging. 
Previous authors have often omitted any form of complete- 



ness correction, instead setting high X-ray flux limits (e.g. 
lUeda et alj|2003h or presenting maximal XLFs assigning all 
unidentified sources to a given redshift bin (e.g. lBarger et aU 
2005). Silv erman et alj (|2008l ) did apply a completeness cor- 
rection, based on the fraction of sources with a given X-ray 
flux with redshift information, although samples were re- 
stricted to objects with bright optical c ounterparts (r ' < 24 ) 
which may introduce additional biases. lEbrero et alj (|2009l ) 




Figure 9. XLF results for the LADE model, constrained at high redshifts using the colour pre-selected samples as described in section 
13.3.21 Binned estimates at z < 1.2 use the hard-band selected sample only, including photometric redshift information (black circles). In 
the final 2 panels binned estimates based on the high-redshift colour pre-selected samples that are used to constrain the evolution are 
shown (red triangles) . At z > 1.2 we also show binned estimates using the hard-band selected sample, although these data were not used 
to fit the model (open circles). The LADE model evaluated at z = is shown in each panel (black dashed lines). 



restricted their analysis to samples with high-redshift iden- 
tification completeness, and applied small additional correc- 
tions based on the fraction of identified sources as a function 
of X-ray flux, although the completeness was likely to be 
over estimated in the CDF-S as they did not consider the 
high fraction of spurious coun t erpart s. A similar approach 
was adopted by lYencho et alj (|2009l ) . although their sam- 
ples were restricted to spectroscopic identifications and thus 



these completeness corrections were large (~ 60%) at the 
faintest flux levels, and may be subject to additional biases 
and inaccuracies given the large number of complex factors 
involved in spectroscopic success rates. We thus believe our 
approach of considering only robust optical counterparts to 
X-ray sources and modelling the distribution of fx/ f opt ra- 
tios, in addition to our colour pre-selection approach at high- 
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redshifts, constitutes a more robust treatment of complete- 
ness than in previous work. 

In Figure [10] we directly compare the evolution of the 
space density of AGN predicted by our LADE model and 
a number of recent studies of the evolution of the XLF, all 
of which concluded that LDDE was taking place. Results 
are shown in two luminosity bins, which approximately cor- 
respond to above and below the characteristic break lumi- 
nosity, Z/„, although L* does evolve and thus such space 
density plots may provide a distorted view of the evolution- 
ary behaviour of the XLF. Indeed, in the higher luminosity 
bin (44.5 < logLx < 46) our model indicates a very dif- 
ferent evolution of the space density at z < 1 to previous 
(LDDE) results, but this is partly due to small offsets in 
luminosity space which manifest as very large differences 
in density because of the steep slope of the bright end of 
the XLF. This can also lead to large differences in binned 
estimates depending on the methodology or assumed evolu- 
tion of the XLF. However, in this regime, the XLF remains 
poorly constrained, and despite the very different evolution- 
ary behaviour our bin ned results are consisten t with the 
ISilverman et al ] <|2008h and lEbrero et a l. (2009) models at 
the ~ 1 — 2a level, although the lYencho et all 1 20091 ) model 
is significantly lower. The high redshift form of our evolu- 
tion in this luminosity bin predicts similar space densities 
to the previous studies, although the form of the evolution 
and location of the peak in the number density is somewhat 
different. 

Our results are in good agreement at low redshifts 
(z < 1) in t he lower luminosity bin (43 < logLx < 44.5), 
although the I Yencho et all (|200gh results remain below our 
estimates. However, at higher redshifts our model predicts 
a muc h higher space density th an found by lYencho et al.l 
(|2009h or ISilverman et all <|2008h . We attribute the decline 
found by these authors to incompleteness, catastrophic fail- 
ures of photometric redshift estimates, and false associations 
of X-ray sourc es with lower redsh ift optical counterparts. 
Conversely, the lEbrero et all l|2009h LDDE model is in very 
good agreement with our result, most likely because of the 
higher completeness of their samples. However, our results 
have shown that an LDDE model is not necessary to describe 
this evolution. 

It is worth noting that our LADE form of the evolu- 
tion is able to describe a "cosmic-downsizing" behaviour in 
which the peak in the space density moves to lower red- 
shifts for lower luminosity AGN, and indeed we do find 
small redshift offsets between the peak space densities for 
different luminosity ranges (e.g. Figure [TO)) . However, our 
model does not find the strongly luminosity-dependent shift 
in peak space density indicated by some p rior studies (e.g. 
lHasinger et al.ll2005l ; [Silverman et alj|200Sl ). This is mainly 
due to our steeper faint-end slope (and thus higher number 
density of low-luminosity AGN at z > 1), but is also due 
to the differences in the form of the bright-end evolution 
discussed above. 



5.2 Remaining uncertainties 

While we believe our work offers a number of improvements 
over previous studies, there are remaining uncertainties as- 
sociated with our approach. 
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Figure 10. The evolution of the space density of AGN based on 
our LADE model for two luminosity ranges: 43 < logLx < 44.5 
(circles) and 44.5 < logLx < 46 (squares). Binned estimates 
are shown from the hard X-ray samples at z < 1.2 (solid sym- 
bols) and the colour pre-selected samples at z ~ 2-3 (open sym- 
bols). The space de nsity based on L DDE models determ ined by 
lEbrero et al] ll2009l. dot-dashed Z me), lYencho et al.l (f2009l . dashed 
line) and lSilverman et al] ]2008l . dotted line) are also shown. 



5.2.1 Photometric redshifts 

As previously discussed, the introduction of photometric 
redshifts, essential to improve the completeness of samples 
at faint X-ray fluxes, does introduce significant uncertainty, 
which we have tried to account for by considering probability 
distributions for the redshift estimates. However, there are 
a number of limitations to photometric redshift techniques 
that are thus also limitations in our studies. This is partic- 
ularly apparent in the catastrophic failure of our photo-z 
estimates at z > 1.2, which we believe led to a biased result 
for the evolution of the XLF, favouring an LDDE scheme 
with a flattening faint-end slope (section l4.2l tabled . Addi- 
tionally, the majority of our photo- z estimates are based on 
template fitting approaches that utilise a variety of galaxy 
templates (and priors based on the galaxy population), and 
thus the derived redshift probability distributions may not 
correspond to the true redshift for an AGN, which could bias 
our results at low as well as high redshifts. Incorporating 
AGN templates and using priors based on the AGN popula- 
tion may offer significant improvements and reduce system- 
atic biases and cat astrophic failures (e.g. IZheng et alJl20o3 : 
ISalvato et ahlliooih . Expanding the number and wavelength 
range of the filters, such as including near- and far-IR data, 
also has potential to improve photometric reds hift estimates 
and p rovide better fits to AGN templates (e.g lSalvato et all 
120081) . 

In addition to the template fitting approach, we also 
utilised ANN2 redshift estimates in the AEGIS field, in 
which an artificial neural network is trained to provide 
photometric redshifts using a sample of objects with spec- 
troscopy. While this approach produced reliable results, it 
is limited to the subset of objects with a certain range of 
colours and magnitudes for which follow-up spectroscopy is 
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available. Indeed, spectroscopic validation of photo-z esti- 
mates for large numbers of optically-faint X-ray sources is 
impossible with current instrumentation, which is a major 
limitation of photometric redshift techniques. 

5.2.2 High-z colour pre- selection 

To reduce the bias at high redshifts we incorporated rest- 
frame UV colour pre-selected samples. This should avoid the 
systematic biases associated with catastrophic failure of the 
photometric redshift estimates, but does risk introducing ad- 
ditional uncertainty and biases. As noted in section l4~2l the 
binned estimates based on the colour pre-selected samples 
at 2.0 < z < 2.5 (dominated by BX selected objects) fell sys- 
tematically below the model, ostensibly indicating a rise in 
the space density of AGN between z ~ 2 and z ~ 3 (also seen 
in Figure [10)) . The BX criteria select objects in a very thin 
slice of colour space, identifying objects with a specific spec- 
tral shape rather than strong spectral features. This makes 
the BX selection both more susceptible to contamination 
by lower redshift sources and more prone to incompleteness 
than the more widely-used LBG selection that identifies a 
strong break in the spectrum shortward of the Lyman-limit. 
The shape and efficiency of the derived selection functions 
will be highly dependent on the modelling of the AGN spec- 
tra and the scatter in their properties, and it is quite possible 
that our relatively simple modelling may have led to over- 
estimates of the completeness. Indeed, both the LBG and 
BX selections require AGN to be bright in the rest-frame 
UV, and will therefore miss low- luminosity, moderately or 
hea vily obscured AG N residing in red, evolved host galax- 
ies. lAird et all (|2008l ) found that LBG selection recovered 
a number density of low-luminosity AGN that is at least as 
high as found via direct X-ray selection and follow-up, tenta- 
tively indicating that the AGN population at z ~ 3 is dom- 
inate d by objects with blue colours (cf. z ~ l. lNandra et alj 
120071 '). However, this may not be true at slightly lower red- 
shifts, resulting in a lower number density at 2.0 < z < 2.5 
via the BX selection as red galaxies will be missed. 

As discussed in section 12.2.41 our combination of LBG 
and X-ray selection to identify z ~ 3 AGN does not suffer 
from low-redshift contamination, but this is not the case for 
the BX selection. Without complete spectroscopic follow-up 
of the X-ray detected BX candidates we cannot be confident 
that we have excluded all interlopers (although we believe 
we are able to identify low-redshift star-forming galaxies, see 
Figure [4|, and it is unclear how this might affect our derived 
evolution of the XLF. Extensions of the BX selection tech- 
nique, which incorporate data at redder wavelengths, may 
allow low-redshift interlopers to be excluded, and help iden- 
tify AGN in redder host galaxies. Indeed, a hybrid of the 
photometric redshift and colour pre-selection approaches, 
utilising data from all available wavebands to perform SED 
fitting and estimate a redshift, but also performing a num- 
ber of colour cuts to avoid catastrophic failures, and care- 
fully modelling the distribution of AGN spectral properties 
to determine the incompleteness, may be the optimal ap- 
proach to determine the high-redshift evolution. 

It is also worth noting that our approach of combin- 
ing direct hard X-ray selection at z < 1.2 and colour pre- 
selection at high redshifts provides no data for z ~ 1.2 — 2.0. 
Significant evolution may be taking place over this redshift 



range , particularly at high luminosities (e.g. lSilverman et al.l 
l2005f ). Indeed, at high luminosities incompleteness issues 
should be less severe, although we note that at z ~ 2 we 
expect ~10% of AGN with L x = 10 445 ergs" 1 will have 
optical magnitudes R > 24, based on our fx/ ) opt distribu- 
tions. In this redshift range, photometric redshift estimates 
are prone to catastrophic failure, and spectroscopic identifi- 
cation can also be difficult even for bright targets. Thus care 
must be taken to avoid biased results. As found in section 
14.21 such systematic biases can erroneously indicate the need 
for a more complex evolutionary model. 

5.2.3 Completeness corrections 

To account for objects with optical magnitudes below our 
limit of R = 25.5, we introduced a simple completeness 
correction assuming the log/x// op t distribution can be de- 
scribed by the same Gaussian function at all redshifts (sec- 
tion [3]2}. This may be an over simplification, as the red- 
shift distribution of optically-f aint sources wi l l diffe r from 
the brighter population. Indeed TlVIainieri et al.l l |200 5l found 
that optically-faint X-ray so urces tended to be foun d at 
higher redshifts. Additionally, iKoekemoer et alj <|2004h dis- 
cuss how Extreme X-ray Objects (EXOs) which are X-ray 
bright but lack optical counterparts, even in very deep Hub- 
ble ACS imaging, may represent a population of very high- 
redshift [z ~ 6 — 7) AGN where the redshifted Lyman- break 
suppresses all emission at optical wavelengths. More sophis- 
ticated modelling of the rest-frame UV-optical properties of 
the AGN population and the effects on observed magnitudes 
as a function of redshift would enable improvements to our 
completeness corrections. Furthermore, the shape and evo- 
lution of the optical luminosity function of the AGN host 
galaxies will affect the completeness of our samples , could 
alter the shape and efficiency of the optical selection func- 
tions, and may systematically impact the photo- z estimates. 
In addition, we do not currently utilise any information from 
the optically-unidentified population (their X-ray fluxes, the 
fact that they lack counterparts to the limits of our imag- 
ing), which could place further constraints on the evolution 
of the XLF. 

5.2.4 Cosmic variance 

Our investigation of the XLF uses an unprecedented large 
area of deep 200 ks Chandra data in the AEGIS field, which 
supplements the data f rom the Chandra Deep Fields used 
in pr evious studies (e.g. iBarger et al. 2005; Silve rman et al.l 

2008) . As such, our determination of the faint-end evolu- 
tion should be less susceptible to the effects of cosmic vari- 
ance than in prior work. However, cosmic variance could in- 
troduce additional uncertainty in our XLF determinations 
(and thus weaken evidence for more complex evolutionary 
schemes). Indeed, assuming our X-ray detected AGN clus- 
ter according to a power-law with scale length ro = 7.1 Mpc 
and s lope 7 = 1.8 (typical of X-ray AGN at z ~ l lCoil et al.l 

2009) , we predict a fractional RMS variation of ~ 5 — 10% 
in the number density of AGN. This is comparable to our 
Poissonian uncertainties due to the number of objects in 
bin, although the errors will be basically covariant across 
the luminosity bins, and thus will not affect the shape of 
the XLF. 
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5.2.5 Intrinsic absorption 

We have calculated rest-frame 2-10 keV luminosities for 
our sources from their hard band fluxes (or soft band for 
the high-z sample), assuming a power-law with T = 1.9. 
This approach should correct for absorption effects due to 
intrinsic column densities of Nh < 10 23 cm -2 at z ~ 1. 
However, the effectiveness of this correction varies with red- 
shift, and this may bias our measurements of the luminos- 
ity, and alter the effective sensitivity of our observations. 
The extent of these effects will depend on the distribution 
of AGN absorption properties, which may be evolving with 
redshift or have a lum inosity dependence l|Ueda et alj 2003; 
La Franca et al.ll2005l). Some previous studies ( Ueda et "all 
20031 ; lLa Franca et al.ll2005l ; lEbrero et all [20091 ) instead at- 
tempted to correct for the effects of absorption on a source- 
by-source basis, and directly account for the Nh distribution 
in their calculation of the XLF. Determining the extent of 
absorption is difficult for our faint sources, which have few 
counts and thus little X-ray spectral information. The effects 
of absorption on our X-ray sensitivity and completeness cal- 
culations will be complicated, and should be corrected for a 
number of potential selection biases . Such effects are in fact 
likely to introduce considerable additional uncertainty as to 
the evolution of the faint end of the XLF, and could thus 
strengthen our assertions of the lack of a requirement for a 
more complicated LDDE model to describe the evolution of 
the XLF. We defer these considerations to future work. 
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5.3 The evolution of AGN accretion activity 

Our determination of the evolution of the XLF presented in 
this paper sheds new light on how the distribution of AGN 
accretion activity evolves over the history of the Universe. 
We conclude that given the currently available data, the 
evolution of the XLF can be described by a luminosity and 
density evolutionary scheme, in which the XLF retains the 
same shape, but shifts in both luminosity and density with 
redshift. A more complicated luminosity-dependent scheme 
does not provide a significantly better description of the 
data; even so our best-fit LDDE for the hard X-ray sample 
only, which we expect is biased at high redshifts, indicates 
a steeper faint-end slope and less evolution of the shape of 
the XLF than mos t prev ious work. 

iHopk ins et all (|2005f ) presented an interpretation of the 
AGN luminosity function in which the bright end traces the 
mass distribution of black-holes accreting at the peak rate 
over their lifetimes, while the faint end corresponds to AGN 
during transitionary periods as they appro ach or decline 
from t heir peak rate of activity. Furthermore, IHopk ins et al.l 
(2006) proposed that a luminosity-dependent lifetime for 
AGN could explain any observed flattening of the faint- 
end slope, as less time is spent in the transitionary pe- 
riod by the more luminous AGN found at high redshifts. 
Our results, however, show that the shape of the XLF does 
not change, thus a strongly luminosity-dependent lifetime 
may not be required. The lack of evolution in the shape 
of the XLF indicates the processes involved in the trigger- 
ing and feeding of AGN remain the same at all redshifts, 
but increase in overall density from the earliest times un- 
til z ~ 1, and moving to lower luminosity systems be- 
tween z « 1 and the present day. Such behaviour could 



Figure 11. 2-10 keV luminosity density of AGN as a function 
of redshift for our LADE model integrated over the luminosity 
ranges indicated. The grey shaded region indicates the la uncer- 
tainty on the total luminosity density of AGN. The luminosity 
density is dominated by moderate luminosity (43 < log L x < 45) 
AGN, and exhibits a mild decline above z as 1.2. 



be consistent with a merger-driven model of the build-up of 
galaxy and black-hole mass, and the triggering of AGN (e.g. 
iKauffmann fc HaehnettfeoOOHHopkins et alj|2008h . in which 
the hierarchical build-up of galaxies increases the number 
density of AGN in the early history of the Universe, but 
the exhaustion of gas supplies in the most massive galaxies 
at later times (e.g. ICowie et al.lll996T ) results in the "down- 
sizing" of AGN to lower luminosities. Fully reconciling such 
a picture with the full range of observational data for AGN, 
in addition to the observed evolution of galaxies, using cur- 
rent theoretical models clearly requires significant further 
investigation. Indeed a number of other processes may play 
a role in triggeri ng and fuelling low lumino sity AGN, such as 
minor mergers (iHernquist fc Mihodll995l ). accretion of gas 
from a hot halo dFa bian & R eeslll995t ) or bar-driven accre- 
tion (|Sellwood fc Moorelll999i r It is essential to assess the 
importance of such processes and understand the role they 
play in the evolution of AGN accretion activity. 

In Figure [11] we plot the 2-10 keV luminosity den- 
sity, f Lx</>(Lx, z)d log Lx, as a function of redshift, which 
provides a tracer of the total AGN accretion activity. The 
evolution of the total luminosity density for all AGN with 
Lx > 10 42 erg s _1 is shown as well as the contribution of 
AGN in set luminosity ranges. The total luminosity den- 
sity peaks at z = 1.2 ± 0.1, rapidly declining to lower red- 
shifts, but the decline to higher redshifts is much milder 
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Figure 12. Total accreted black-hole mass density against look- 
back time based on our LADE {solid, grey shaded region indicates 
la uncertainty in the derived model). The lower curves corre- 
spond to the same luminosity ranges indicated in Figure [TT] We 
find ~ 50% of the local black-hole mass density is built up in 
AGN actively accreting at z > 1. 



indicating significant AGN activity is taking place out to 
z ~ 3-4. AGN with luminosities in the range L x » 10 43 ~ 45 
erg s - (i.e. ~ L») are responsible for the majority of the lu- 



minosity density, although the contribution of Lx = 10 44 45 

1 AGN falls off rapidly at z < 0.8 as L* evolves, 
-V13-44 „_„ „-i 



and 



erg s 

L x = 10 43 - 44 erg s -1 AGN dominate. 

We can track the build-up of black-hole mass more di- 
rectly by rel ating AGN lum inosity to mass accretion, as first 
proposed bv lSoltanl <|l982h . 



eM a , 



eM hh c 2 



(26) 



where M acc is the mass accretion rate, Mbh is the rate 
of change of black-hole mass density, Lboi is the bolomet- 
ric luminosity, e is the radiative efficiency of the accretion 
process and c is the speed of light. We have adopted a 
simple approach (e.g. iBarger et al.l [200ot ). converting our 
2-10 keV X-ray luminosities to bolometric value s using a 
constant conversion factor (40. lElvis et al.lll994u . and as- 
suming a single valu e of the radiative efficiency, e = 0.1 
|Marconi et alj|2004h . We note that a number of authors 
have discussed the ne ed for luminosity-d e pendent bolo- 
metr ic cor rections (e.g. iMarconi et al. 120041 : iHopkins et al.l 



120071 ). and I Vasudevan fc Fabianl ( 2007 ) reported significant 
object to object variation in bolometric corrections, which 
may depend on the Eddington ratio. The radiative efficiency 
may also vary significantly between AGN , depending on the 
spin of the black-hole (e.g. lThornelll974h i n addition to the 
speci fics of the accretion processes (e.g. iMerloni fc Hein5 
2008). However our simple assumptions allow an initial in- 
vestigation of the consequences of our derived XLF evolu- 
tion on the build-up of black-hole mass. Using equation [26] 
we can convert our luminosity densities to a mass accretion 
rate, and thus calculate the total black-hole mass density as 
built up by accretion activity over the history of the Uni- 
verse. Our results are shown in Figure [121 



Based on our LADE model we predict a local black-hole 
mass density of 2.2 ± 0.2 x 10 5 Mq Mpc -3 , where the er- 
ror reflects the uncertainties in our model fit, but not the 
potentially larger uncertainties in bolometric correction or 
accretion effici ency. This value is in good agreement with 
the estimate of lYu fc Tremaind (]2002t ) (2.5 ± 0.4 x 10 5 M & 
Mpc -3 ) based on velocity dispersions of early- type galaxies 
in SDSS and th e Mbh — Q" rela ti on, a lthough is lower than 
the estimate of lMarcon[ et all (|2004l ) x 10 5 M 

Mpc -3 ), possibly indicating that our XLF does not pro- 
vide a complete census of the history of accretion activity. 
Figure [12] shows that a significant fraction (~ 50%) of this 
total mass density is accreted at z > 1. While the major- 
ity of the mass build up takes place in moderate luminosity 
AGN (Lx = 10 43-45 erg s -1 ), a significant fraction is accu- 
mulated at lower luminosities (Lx = 10 42-43 erg s -1 ). The 
LDDE model from section [4] predicts a lower local black-hole 
mass density, mainly due to the smaller numbers of AGN at 
these low luminosities and high redshifts. The redshift range 
z ~ 1-3 clearly corresponds to a period of significant AGN 
activity, and thus it is essential to accurately measure the 
XLF down to Lx ~ iq 42 - 4 3 er g g -i m this epoch to deter- 
mine the history of black-hole mass accretion. 

We can also compare our derived black-hole mass ac- 
cretion rates to star-formation r ates, which we show in 
Figure | 13[ Previous studies (e.g. iBovle fc Terlevichl 1 19981 : 
ISilverman et al . 2008) have shown close similarities between 
the rapid increase in star-formation rate and black-hole ac- 
cretion at z < 1, which we confirm. Our results indicate 
that this correlation continues out to high redshifts (cf. 
ISilverman et al.ll2008h . at l east when comparing to the re- 
cent star- formation rates of lBouwens et alTl|2007l ). although 
we note that our model is extrapolated far beyond the red- 
shift range probed by our data. Comparisons of the galaxy 
and AGN luminosity functions may reveal differences in the 
details of the evolving distributions of activity, which could 
reveal further facets of the co-evolutionary processes and the 
feedback regulating AGN activity and star-formation. 



6 SUMMARY 

We have presented new observational determinations of the 
X-ray luminosity function of AGN. We utilised a hard X- 
ray selected sample from both the 2 Ms Chandra Deep 
Fields and the large area, deep 200 ks Chandra survey in 
the AEGIS field. A likelihood ratio method was employed 
to match X-ray sources to optical counterparts and ensure 
only robust associations were considered. To improve our 
redshift completeness, we supplemented spectroscopic iden- 
tifications with photometric redshifts. However, we found 
our photo- 2 estimates were systematically biased and prone 
to catastrophic failure at z > 1.2. We therefore adopted a 
rest-frame UV colour pre-selection approach to constrain the 
XLF at high redshifts. By carefully modelling the expected 
colour distributions of AGN, performing simulations of the 
optical data and repeating the source recovery procedure, 
we were able to determine well-defined selection functions 
that allowed us to correct for the incompleteness associated 
with the pre-selection technique. 

We developed a sophisticated Bayesian methodology to 
determine the evolution of the XLF, which accounts for the 
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Figure 13. Com parison of star- formation rates (circles) from 
the compilation of [Hopkins (2004|) and black- hole accretion rate 
based on our LADE model (solid line; grey shading indicates la 
uncertainty in the model). Black- hole mass accretion rate assumes 
a constant bolometric correction and accretion efficiency, e = 0.1, 
and is scaled up by a factor 5000 l lSilverman et al.l 120081) . We 
also show rece nt measurements of th e star-formation rate at high 
redshifts from lBouwens et al.l ||2007| . squares), which indicate the 
correlation between star-formation and black-hole accretion rate 
continues at high redshifts. 



uncertainties in photometric redshift estimates, the Poisso- 
nian nature of the X-ray flux estimate, the fraction of sources 
with counterparts below the magnitude limits of the opti- 
cal data, and the optical selection functions at high red- 
shifts. Using a Bayesian model comparison approach, we 
found that a Pure Luminosity Evolution was unable to ad- 
equately describe the evolution of the XLF, but a scheme 
in which the XLF evolves in both luminosity and density 
(LADE), but does not change shape, did provide a good 
fit. We did not find significant evidence for a more com- 
plex luminosity-dependent density evolution (LDDE), and 
the associated flattening of the faint-end slope. However, an 
LDDE model was required if the direct X-ray selected sam- 
ples and the photometric redshift estimates were utilised at 
z > 1.2, but we believe such a result to be biased due to the 
inclusion of unreliable photo- 2 estimates. 

The form of our derived evolu tion of the XLF dif- 
fers from many previous results (e.g. ISilverman et alj 120081 ; 
lEbrero et al.ll2009l ; lYencho et al.ll2009l ). Our XLF retains the 
same shape at all redshifts, which may have a number of con- 
sequences for the evolution of the black-hole mass function, 
AGN lifetimes and duty cycles, and the processes that fuel 
and trigger AGN. We find that the total luminosity den- 
sity peaks at z = 1.2 ± 0.1. While there is a mild decline 
at z > 1.2, significant AGN activity is still taking place at 
high redshifts. We find significantly higher number densities 
of low luminosity AGN (Lx < 10 44 erg s _1 ) at z > 1, which 
make a significant contribution to the to the total luminos- 



ity density and growth of black-hole mass at these redshifts. 
Our results do not indicate a strong shift of the peak space 
density to lower redshifts for lower luminosities, although a 
small shift is predicted by our best-fit evolutionary model. 
The mild decline in AGN activity to high redshifts appears 
to correlate with that seen in the star-formation rate, con- 
sistent with a co-evolutionary scheme for black-holes and 
galaxies. 
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